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Abstract 

We  consider  the  following  maximum  disjoint  paths  problem  (mdpp).  We  are  given 
a  large  network,  and  pairs  of  nodes  that  wish  to  communicate  over  paths  through  the 
network  —  the  goal  is  to  simultaneously  connect  as  many  of  these  pairs  as  possible 
in  such  a  way  that  no  two  communication  paths  share  an  edge  in  the  network.  This 
classical  problem  has  been  brought  into  focus  recently  in  papers  discussing  applications 
to  routing  in  high-speed  networks,  where  the  current  lack  of  understanding  of  the  mdpp 
is  an  obstacle  to  the  design  of  practical  heuristics. 

We  consider  the  class  of  densely  embedded,  nearly- Eulerian  graphs,  which  includes 
the  two-dimensional  mesh  and  many  other  planar  and  locally  planar  interconnection 
networks.  We  obtain  a  constant-factor  approximation  algorithm  for  the  maximum  dis¬ 
joint  paths  problem  for  this  class  of  graphs;  this  improves  on  an  0(logn)-approximation 
for  the  special  case  of  the  two-dimensional  mesh  due  to  Aumann-Rabani  and  the  au¬ 
thors.  For  networks  that  are  not  explicitly  required  to  be  “high-capacity,”  this  is  the 
first  constant-factor  approximation  for  the  MDPP  in  any  class  of  graphs  other  than 
trees. 

We  also  consider  the  MDPP  in  the  on-line  setting,  relevant  to  applications  in  which 
connection  requests  arrive  over  time  and  must  be  processed  immediately.  Here  we 
obtain  an  asymptotically  optimal  0(logn)-competitive  on-line  algorithm  for  the  same 
class  of  graphs;  this  improves  on  an  0(lognloglogn)-competitive  algorithm  for  the 
special  case  of  the  mesh  due  to  Awerbuch,  Gawlick,  Leighton,  and  Rabani. 


1  Introduction 

We  consider  the  following  maximum  disjoint  paths  problem  (mdpp).  We  are  given  a  large 
network,  and  pairs  of  nodes  that  wish  to  communicate  over  paths  through  the  network  - 
the  goal  is  to  simultaneously  connect  as  many  of  these  pairs  as  possible  in  such  a  way  that 
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no  two  communication  paths  share  an  edge  in  the  network.  This  problem  is  well-known  to 
be  computationally  difficult.  Deciding  whether  all  pairs  can  be  so  connected  is  one  of  Karp’s 
original  NP-complete  problems  [12];  it  remains  NP-complete  even  when  the  underlying  graph 
is  the  two-dimensional  mesh  [15]. 

Our  interest  in  this  problem  comes  from  two  main  sources.  First,  establishing  disjoint 
paths  is  fundamental  to  routing  in  high-speed  networks  (see  for  example  the  applications 
mentioned  in  [6,  8,  21],  as  well  as  applications  to  optical  routing  in  [1,  3,  24]).  Although  the 
types  of  routing  problems  that  arise  in  such  settings  tend  to  have  additional  side  constraints 
(e.g.  connections  have  limited  duration  and  can  bring  varying  amounts  of  “profit”),  the 
formulation  described  in  the  first  paragraph  contains  the  essence  of  virtually  all  such  real- 
life  routing  problems  in  which  each  connection  consumes  a  large  fraction  of  the  bandwidth 
on  a  link.  As  such,  the  current  lack  of  understanding  of  the  disjoint  paths  problem  is  a 
major  obstacle  to  the  design  of  practical  heuristics.  Indeed,  [6]  notes  that  in  practice,  the 
greedy  algorithm  tends  to  be  used  for  routing,  despite  its  bad  performance  on  a  number 
of  very  common  interconnection  patterns.  Moreover,  robust  ways  are  known  for  converting 
algorithms  for  the  MDPP  into  algorithms  that  can  handle  connections  of  limited  duration  or 
variable  value  [5];  thus,  the  difficulties  contained  in  these  more  elaborate  routing  problems 
seem  to  stem  mainly  from  the  intractability  of  the  MDPP. 

This  problem  is  also  of  basic  interest  in  algorithmic  graph  theory.  A  lot  of  work  has 
been  done  on  identifying  special  cases  of  the  disjoint  paths  problem  that  can  be  solved  in 
polynomial  time,  or  for  which  simple  min-max  conditions  can  be  stated;  see  the  survey  by 
Frank  [10].  Much  less  work  has  been  done,  however,  on  approximation  algorithms  for  the 
MDPP;  we  are  interested  in  extending  the  classes  of  graphs  for  which  good  approximations 
can  be  obtained. 

1.1  Our  Results 

To  be  precise,  let  G  =  (V,  E)  be  a  graph  on  n  vertices  and  T  =  {(si,  ti), . . . ,  (s^,  f*,)}  a 
collection  of  terminal  pairs  —  pairs  of  vertices  of  G.  We  say  that  T  is  realizable  in  G  if 
there  exist  mutually  edge-disjoint  Sj-R  paths,  for  i  —  1, . . . ,  k.  The  problem  is  then  to  find  a 
realizable  subset  of  T  of  maximum  cardinality. 

Our  first  main  result  is  a  constant-factor  approximation  for  the  maximum  disjoint  paths 
problem  in  the  class  of  densely  embedded,  nearly- Eulerian  graphs  (defined  below),  which 
includes  many  common  planar  and  locally  planar  interconnection  networks.  This  improves 
on  an  O  (log  ^-approximation  for  the  case  of  the  two-dimensional  mesh  due  to  Aumann  and 
Rabani  [3]  and  an  O  (log  ^-approximation  for  a  class  of  planar  graphs  including  the  mesh 
due  to  the  authors  [14],  Our  present  algorithm  makes  use  of  variants  of  a  number  of  the 
techniques  developed  in  our  earlier  paper  [14], 

The  assumption  that  we  know  all  the  terminal  pairs  in  advance  is  not  reasonable  in 
situations  in  which  connection  requests  between  pairs  of  nodes  arrive  over  time  and  must 
be  processed  immediately.  In  such  a  setting,  it  makes  sense  to  consider  on-line  routing 
algorithms.  Such  an  algorithm  is  given  the  graph  G,  terminal  pairs  arrive  in  an  arbitrary 
order,  and  for  each  such  pair  it  must  irrevocably  reject  it,  or  assign  it  a  path  in  G.  As 
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is  standard,  we  refer  to  the  approximation  ratio  achieved  by  an  on-line  algorithm  as  its 
competitive  ratio ;  such  an  algorithm  is  said  to  be  c-competitive  if  its  competitive  ratio  is  at 
most  c. 

Our  second  main  result,  then,  is  an  O (log  n)-competitive  randomized  on-line  algorithm  for 
the  mdpp  in  densely  embedded,  nearly- Eulerian  graphs.  This  improves  on  an  0(log  n  log  log  n)- 
competitive  algorithm  for  the  special  case  of  the  two-dimensional  mesh  due  to  Awerbuch, 
Gawlick,  Leighton,  and  Rabani  [6];  moreover,  [6]  proves  that  no  randomized  on-line  algo¬ 
rithm  for  the  two-dimensional  mesh  can  be  better  than  Q(log  n)-competitive,  implying  that 
our  algorithm  is  asymptotically  optimal. 

We  feel  that  an  important  feature  of  our  algorithms,  in  addition  to  the  improved  bounds, 
is  that  they  are  not  specific  to  the  mesh;  the  advantage  of  developing  algorithms  that  work 
on  the  somewhat  larger  class  of  “densely  embedded  graphs”  is  that  are  not  sensitive  to  small 
variations  in  the  structure  of  the  underlying  graph.  This  could  be  of  value  in  the  context  of 
network  routing,  where  the  underlying  network  may  have  a  “mesh-like”  topology,  but  lack 
the  completely  regular  structure  of  the  mesh.  In  contrast,  previous  algorithms  such  as  [6,  3] 
could  not  be  applied  to  any  network  other  than  the  mesh  itself,  since  they  required  its  fixed 
row/column  structure. 

The  size  of  the  constants  in  our  algorithms  as  presented  here,  while  not  astronomical, 
pushes  them  outside  the  range  of  immediate  practical  utility.  However,  the  previous  best 
bounds  —  both  off-line  and  on-line  —  for  the  two-dimensional  mesh  [3,  6,  14]  involve  similarly 
large  constants  inside  the  ()(■)  notation.  Moreover,  despite  the  large  constants,  some  of  the 
ideas  used  by  the  algorithms  here  may  be  of  use  in  suggesting  practical  heuristics. 

The  rest  of  the  paper  is  organized  as  follows.  First,  we  present  some  preliminary  algo¬ 
rithmic  tools  that  we  need  in  Section  2.  Our  routing  algorithms  are  simpler  to  explain  in 
the  special  case  of  the  mesh,  and  we  consequently  present  this  case  first.  This  will  be  done 
in  Section  3.  In  Section  4,  we  then  introduce  the  class  of  densely  embedded  graphs,  and 
we  present  the  algorithms  for  this  class  in  Section  5.  Finally,  we  show  how  to  extend  our 
algorithms  to  a  slightly  more  general  class  of  graphs  in  Section  6. 

1.2  Previous  Work 

Much  of  the  previous  work  on  this  problem  has  dealt  with  the  case  in  which  each  path 
consumes  only  a  small  fraction  of  the  available  bandwidth  on  an  edge;  this  can  be  modeled 
by  requiring  fi(logn)  parallel  copies  of  each  edge.  In  this  case,  the  randomized  rounding 
technique  of  Raghavan  and  Thompson  [23,  22]  can  be  used  to  obtain  an  off-line  constant- 
factor  approximation.  Awerbuch,  Azar,  and  Plotkin  give  an  on-line  O  (log  n)-competitive 
algorithm  for  this  case  [4],  which  they  show  is  asymptotically  tight  for  deterministic  on-line 
algorithms. 

As  noted  in  [6]  however,  there  are  many  applications  in  which  each  communication  path 
consumes  a  large  fraction  of  the  available  bandwidth  on  a  link;  thus  it  makes  sense  to  con¬ 
sider  approximation  algorithms  for  graphs  without  a  large  number  of  parallel  edges.  The 
results  here  are  much  more  restricted.  For  trees  with  parallel  edges,  Garg,  Vazirani,  and  Yan- 
nakakis  [11]  obtain  an  off-line  2-approximation  (the  maximization  problem  is  NP-complete, 
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though  deciding  realizability  is  easy);  Awerbuch  et.  al.  [6]  give  an  O (log d)-competitive  ran¬ 
domized  on-line  algorithm  for  trees  of  diameter  d,  extending  an  earlier  result  of  Awerbuch 
et.  al.  [5].  Essentially  the  only  approximation  results  known  for  graphs  other  than  trees  are 
those  mentioned  earlier  for  the  mesh  and  related  planar  graphs  [3,  6,  14].  Thus  our  result 
here  is  the  first  constant-factor  approximation  for  any  class  of  graphs  other  than  trees,  when 
one  does  not  require  fl(logn)  parallel  copies  of  each  edge. 

A  different  approach  is  taken  in  papers  of  Peleg  and  Upfal  [21]  and  Broder,  Frieze,  and 
Upfal  [8]  (see  also  Broder  et.  al.  [7]).  Here  the  underlying  graph  G  is  assumed  to  have 
strong  expansion  properties;  in  this  case  one  can  prove  that  any  set  of  terminal  pairs  of  at 
most  a  given  size  must  be  realizable  in  G,  and  that  corresponding  paths  can  be  found  in 
(randomized)  polynomial  time.  The  results  in  [8]  are  strong  enough  that  they  implicitly 
provide  a  polylogarithmic  approximation  for  the  MDPP  in  sufficiently  strong  expanders  of 
bounded  degree. 

In  this  context,  it  is  worth  mentioning  the  following  closely  related  routing  problem: 
one  must  route  all  terminal  pairs  so  as  to  minimize  the  maximum  congestion  on  any  edge; 
that  is,  the  maximum  number  of  paths  that  contain  an  edge.  Deciding  whether  T  can  be 
routed  with  congestion  equal  to  1  is  the  same  as  deciding  whether  T  is  realizable;  but  as  an 
optimization  problem,  minimizing  congestion  is  much  better  understood.  Aspnes  et  al.  [2] 
give  an  on-line  O  ( log  n)-competit.ive  algorithm  for  the  problem;  and  a  randomized  rounding 
algorithm  of  Raghavan  and  Thompson  [23]  gives  a  routing  off-line  with  congestion  at  most 
OPT  +  o(OPT)  +  O(logn),  where  OPT  denotes  the  optimum  congestion.  This  leaves  open 
the  question  of  whether  a  constant-factor  approximation  is  achievable  also  for  small  values  of 
OPT.  In  a  companion  paper  with  Satish  Rao  [13],  we  give  a  constant-factor  approximation 
for  densely  embedded  graphs;  it  was  here  that  the  idea  of  using  a  “simulated  network”  (see 
Section  3)  was  initially  developed. 

Cases  in  which  the  MDPP  can  be  solved  in  polynomial  time  are  surveyed  in  [10];  here 
we  only  discuss  two  specific  results  that  we  will  use  in  handling  densely  embedded  graphs. 
First,  suppose  G  is  planar,  the  terminals  T  lie  on  a  single  face  of  G,  and  the  pair  (G,  T) 
satisfies  the  following  parity  condition :  the  augmented  graph  formed  by  adding  to  G  the  edges 
corresponding  to  T  must  be  Eulcrian.  In  this  case,  a  theorem  of  Okamura  and  Seymour  [20] 
says  that  the  realizability  of  T  in  G  can  be  decided  in  polynomial  time;  and  in  fact  the 
following  cut  condition  is  sufficient  for  realizability:  one  cannot  remove  j  edges  from  G  and 
separate  more  than  j  terminal  pairs.  A  linear-time  algorithm  for  this  problem  has  recently 
been  obtained  by  Wagner  and  Weihe  [30].  We  will  use  an  extension  of  the  Okamura-Seymour, 
due  to  Frank  [9],  which  concerns  the  case  in  which  the  parity  condition  does  need  not  to 
hold  on  the  face  containing  the  terminals. 

We  also  use  a  theorem  of  Schrijver  [28]  that  provides  an  algorithm  for  finding  vertex- 
disjoint  paths  in  a  graph  embedded  on  a  compact  surface  E,  such  that  the  paths  satisfy  given 
homotopy  constraints. 
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2  Preliminary  Tools 

2.1  The  AAP  Algorithm 

We  make  use  of  a  variant  of  an  on-line  mdpp  algorithm  of  Awerbuch,  Azar,  and  Plotkin  [4],  If 
H  is  a  graph  with  n  nodes  in  which  each  edge  has  capacity  at  least  log  2n,  the  algorithm  of  [4] 
achieves  a  competitive  ratio  of  2  log  4n.  For  onr  purposes,  we  need  to  develop  a  strengthening 
of  this  “AAP  algorithm”:  we  want  to  be  competitive  against  the  fractional  optimum;  and 
when  we  deal  with  the  more  general  case  of  densely  embedded  graphs,  we  want  only  to 
require  capacities  to  be  £  log  n,  for  an  arbitrary  £  >  0.  Here,  we  show  how  to  obtain  such  an 
algorithm. 

Proposition  2.1  If  all  edge  capacities  are  at  least  (s  log  n  +  1  +  £),  there  is  a  deterministic 
on-line  mdpp  algorithm  that  is  0(21/£  log n)- competitive  against  the  fractional  optimum. 

Proof.  We  follow  the  AAP  algorithm  and  its  analysis  very  closely.  We  vary  a  little  from 
their  notation,  since  we  only  deal  here  with  routing  a  maximal  number  of  requests,  each  of 
infinite  duration.  Thus,  the  ith  request  is  specified  by  a  pair  of  terminals.  We  define 

the  “profit”  of  the  connection  to  be  n;  thus  the  total  profit  obtained  by  the  on-line  algorithm 
is  simply  n  times  the  number  of  terminal  pairs  routed. 

Define  p  =  2 1+1/£n,  so  we  have 

£  log  p  =  £  log  n  +  1  +  £. 

Let  ue  denote  the  capacity  of  edge  e;  thus  we  can  assume  that  for  all  e, 

ue  >  £  log  p. 

With  this  value  of  p,  we  now  run  the  AAP  algorithm  —  for  the  sake  of  completeness,  we 
state  this  algorithm  here. 

For  j  =  1,2,...,  k: 

Define  \Je  to  be  the  fraction  of  ue  consumed  by  paths  already  routed. 

Define  cfe  =  ue(/iA=  —  1). 

For  a  request  route  it  on  any  path  P  satisfying  J2e&p  y~ce  —  n- 

If  no  such  path  is  available,  then  reject  the  request. 


First  we  argue  why  the  relative  load  on  an  edge  will  never  exceed  1.  At  the  moment 
before  this  happened,  on  edge  e  say,  we  had 


\i  >  i  -  —  >  i  — 

Ue 


1 

£  log  p 
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Thus 


c? 

Ue 


=  fl 


.A  i 


>  fll  elogM  —  1 


2l/e 
>  n. 


1  =  2n  -  1 


So  we  have 

Cg 

—  >  n 
ue 

and  thus  the  connection  could  not  have  used  this  edge. 

Suppose  there  are  a  total  of  k  requests.  Let  A  denote  the  set  of  requests  routed  by  the 
AAP  algorithm.  Then  we  show 


2'+'/- log  (!) 

3  SA  e 

As  in  the  proof  in  [4]  we  show  this  by  induction  on  the  number  of  admitted  requests,  via 
proving  that  if  the  algorithm  admits  the  jth  request,  we  have 

y  4+1  ~  cl  <  2 1+1/£n  log  n . 

e 

So  consider  edge  e  on  the  jth  path  used  by  the  AAP  algorithm.  We  have 

4+1  -4  =  ue  (/tA»( 2(log^/“e  -  1))  . 

Now  the  exponent  on  the  2  is  at  most  1/e,  and  for  x  €  [0, 1/e]  we  clearly  have  2X  —  1  <  21/£-x. 

Thus 


ci+1  -  4  <  ue  ■  /iAe  •  21/£  •  (log  fi)/ue 


=  tixi  •  21/£  ■  log/i 
=  21/,£  •  log  /i  • 


c3 

—  +  1 

ue 


Summing  over  all  edges  gives  the  desired  bound. 
Finally,  we  show  that  the  expression 


E  cj+1  (2) 

e 

is  an  upper  bound  on  the  profit  of  the  fractional  optimum  minus  the  profit  of  the  on-line 
algorithm.  ([4]  shows  this  for  the  integer  optimum,  but  the  proof  is  essentially  the  same.) 
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Let  Q  denote  the  set  of  indices  which  were  rejected  by  the  on-line  algorithm  but  for 
which  a  positive  fraction  of  the  demand  was  routed  by  the  optimum.  For  j  G  Q,  suppose 
that  the  fractional  optimum  uses  paths  Pj , . . . ,  P?,  with  weights  yj , . . . ,  7?.  Then  since  j 
was  rejected  by  the  on-line  algorithm,  and  the  costs  are  monotonic  in  the  indices,  we  must 
have 


n 


<E 


_ 

Ue 


for  any  i,j.  Then  for  any  edge  e  we  have 


E  . 

i,j:eeP '] 


71 

Up 


<1, 


and  hence  we  have 


EEt>  <  EEE 


V-cfc+1 


j  * 


J  *  eSP* 


■uP 


<  E 

e 

<  E 


^+1 


fc+i 


E 
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Combining  the  bounds  in  Equations  (1)  and  (2),  we  obtain  the  claimed  competitive  ratio. 


A  lower  bound  of  [4]  implies  that  the  factor  of  21,/e  is  unavoidable  for  deterministic  on-line 
algorithms. 

2.2  Combining  On-Line  Algorithms 

In  routing  connections  on-line,  we  will  adopt  an  approach  in  which  the  decision  whether  to 
accept  a  given  connection  is  made  by  a  combination  of  several  algorithms  —  the  connection 
is  accepted  if  each  of  the  individual  algorithms  accepts  it.  From  the  competitive  ratios  of 
these  individual  algorithms  one  can  infer  a  competitive  ratio  for  this  combined  algorithm;  in 
this  section  we  show  how  this  can  be  done. 

Let  U  denote  a  finite  set,  with  Si, . . . ,  Sn  subsets  of  U  such  that  U  =  UjA*.  Let  Ti  denote 
a  collection  of  subsets  of  S{  closed  with  respect  to  inclusion,  and  let 

T  =  {C:  Vi  {CnSi)  6T,}. 

Given  a  set  U'  C  U,  define  p(G/)  to  be  the  maximum  size  of  a  member  of  T  contained  in  U' . 
We  wish  to  design  an  algorithm  for  the  following  on-line  maximization  problem  with  respect 
to  U  and  T.  Elements  of  some  U'  C  U  arrive  in  an  arbitrary  order,  and  on  each  element  our 
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algorithm  either  accepts  or  rejects  it;  the  goal  is  to  accept  a  subset  of  these  elements  that  is 
in  T  and  as  large  as  possible  relative  to  Our  algorithm  will  be  called  c-competitive  if 

it  always  accepts  a  set  of  size  at  least  ^n(U'). 

We  can  define  the  corresponding  on-line  maximization  problems  with  respect  to  St  and 
Ti,  for  each  i  =  1 , ,n,  in  exactly  the  same  way.  Say  that  for  each  i,  we  are  given  an 
algorithm  A,  which  is  Q-competitive  for  the  problem  associated  with  St  and  T%.  Moreover, 
we  assume  that  the  state  of  Aj  is  completely  determined  by  the  set  of  elements  it  has  accepted 
so  far.  We  then  define  our  “combined  algorithm”  A  =  A"=  l  At  for  the  on-line  maximization 
problem  with  respect  to  U  and  T  as  follows.  As  each  u  G  U'  is  presented  to  A,  it  accepts  u 
iff  for  each  i  such  that  u  G  Sj,  Aj  accepts  u.  The  total  set  accepted  so  far,  intersected  with 
Si,  serves  as  the  state  for  A*.  Let  c*  denote  the  maximum  competitive  ratio  of  any  of  the 
algorithms  A,,  and  suppose  each  element  of  U  appears  in  at  most  d  of  the  St. 

Proposition  2.2  A  is  c*d- competitive. 

Proof.  Assume  the  algorithm  A  was  presented  with  a  set  U'  and  it  returned  X.  Let  Y  denote 
a  member  of  T  contained  in  U'  of  maximum  size;  we  show  that  |Y|  <  c*d\X\.  Let  R[  denote 
the  elements  of  Y  \  X  that  were  rejected  by  algorithm  Aj,  Jj  =  X  D  Y  fl  Si  the  elements  of 
Si  accepted  by  both  A  and  the  optimal  solution  Y,  and  Ri  =  U  Rf  Note  that  Y  =  U iRif 
and  R,  C  Y  fl  Si  G  T%. 

We  want  to  prove  that  \Ri\  <  c*\X  fl  Si\  (i  —  1, . . . ,  n).  Set  U-  =  (X  P\  Si)  U  Rp  these  are 
the  elements  of  U'  fl  Si  either  accepted  by  A  or  rejected  by  Aj.  Order  U[  as  it  appears  in  U' , 
and  present  it  as  input  to  Aj.  Then  as  in  the  running  of  the  combined  algorithm  A,  Aj  will 
accept  precisely  the  set  X  fl  Si.  Since  A,:  is  c*-competitive,  and  Ri  G  we  have 


\Ri\  <  C*\x  n  Si\. 

We  also  have  |Y|  <  and  by  the  definition  of  d  we  have  \X  fl  Sf  <  d\X\.  Thus 

\Y\  <  E  1^1  <  c*  E  lX  n  <  c*d\x\.u 

i  i 


In  the  natural  way,  one  can  define  a  fractional  version  of  this  model;  one  then  shows  the 
following  by  a  minor  modification  of  the  proof  above. 

Proposition  2.3  If  each  Aj  is  c* -competitive  against  the  fractional  optimum,  then  A  is 
c* d- competitive  against  the  fractional  optimum. 

3  The  Two-Dimensional  Mesh 

In  this  section,  let  G  =  ( V ,  E)  denote  the  n  x  n  mesh.  A  very  rough  sketch  of  the  algorithms 
is  as  follows.  Since  much  stronger  results  are  known  for  cases  in  which  edges  have  capacity 
fi(logn),  we  want  to  model  G  by  a  high-capacity  “simulated  network”  J\f .  To  do  this  we 
choose,  for  a  constant  7,  a  maximal  set  of  7  log  n  x  7  log  n  subsqures  of  G  subject  to  the 


condition  that  the  distance  between  subsquares  is  at  least  27  log  n.  These  will  serve  as  the 
nodes  of  A/".  We  now  label  some  pairs  of  subsquares  as  “neighbors”  and  connect  them  with 
fi(log  n)  parallel  edges;  these  are  the  high-capacity  edges  of  A f . 

On  this  network  J\f  we  run  the  algorithm  of  Awerbuch,  Azar,  and  Plotkin  [4]  in  the 
on-line  case,  and  the  algorithm  of  Raghavan  and  Thompson  [23]  in  the  off-line  case.  The  al¬ 
gorithms  running  in  A/"  will  return  routes  consisting  of  a  sequence  of  neighboring  subsquares. 
To  convert  such  a  route  into  a  path  in  G,  we  construct  disjoint  paths  between  neighboring 
subsquares.  We  link  a  sequence  of  neighboring  pairs  together  by  the  natural  crossbar  struc¬ 
tures  surrounding  each  subsquare.  This  leaves  us  with  the  problem  of  routing  from  each 
original  terminal  to  the  boundary  of  its  subsquare.  In  the  on-line  case  we  will  route  at  most 
one  terminal  in  each  subsquare,  and  hence  routing  out  of  a  subsquare  will  be  easy.  In  the 
off-line  case  we  use  network  flow  techniques  to  route  the  appropriate  subset  of  terminals  to 
the  boundary  of  the  subsquare.  Finally,  to  prove  the  approximation  ratios,  we  argue  that 
the  number  of  pairs  routed  by  the  optimum  in  G  is  upper-bounded  by  the  maximum  number 
of  pairs  that  can  be  routed  in  a  copy  of  J\f  in  which  all  edges  have  capacity  O  ( log  n) . 

When  G  denotes  the  two-dimensional  mesh,  let  G[i,j]  denote  the  vertex  with  row  number 
i  and  column  number  j,  and  G[i  :  i',j  :  j']  the  subsquare 

{G\p,q\  ■ i<p<i ',  j  <q<j’}. 

Let  d(u,v )  denote  the  least  number  of  edges  in  a  u-v  path.  By  the  L distance  between 
vertices  G[i,j]  and  G[i',j'],  we  mean  L00(G[i,j],G[i',jr\)  =  max(|i  —  i! |,  | j  —  j'\). 

3.1  Building  the  Simulated  Network 

We  choose  a  constant  7  >  1  (any  constant  will  do;  it  will  have  an  influence  on  the  approx¬ 
imation  ratio  we  obtain).  Our  first  goal  is  to  choose  a  maximal  set  of  “mutually  distant” 
vertices  around  which  to  grow  nodes  of  the  simulated  network.  We  divide  the  the  mesh  into 
7  log  n  by  7  log  n  subsquares  as  follows. 

Definition  3.1  A  subsquare  of  V  is  called  a  7-block  if  for  some  natural  numbers  i  and  j , 
it  is  equal  to  the  set 


G[(i  —  1)7 log n  :  77  log ri,  (j  —  1 ) 7 log 71  :  j'7  log  n] . 

If  X  is  a  7 -block  with  associated  natural  numbers  i  and  j ,  then  the  vertex 

1  1 

V  =  G[{i  -  -blog 77,  ( j  -  ~b log 77] 

will  be  called  the  center  of  the  7 -block,  and  we  will  denote  X  by  Cv.  A  boundary  vertex  of  X 
is  one  with  maximal  or  minimal  row  or  column  number.  We  use  X*  to  denote  the  union  of 
X  with  the  (at  most)  eight  other  'y -blocks  that  share  boundary  vertices  with  X. 
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Figure  1:  A  cluster  and  its  surroundings 

By  a  wall  of  X ,  we  mean  a  maximal  set  of  boundary  vertices  having  the  same  row  or 
column  number.  A  vertex  of  X  is  internal  if  it  is  not  a  boundary  vertex. 

Let  V'  denote  the  set  of  all  centers  of  7-blocks.  We  build  a  graph  G'  on  V'  by  joining 
u,  v  G  V'  if  the  corresponding  sets  C*  and  C*  intersect  at  an  internal  vertex.  We  now  run  a 
randomized  version  of  Lnby’s  maximal  independent  set  algorithm  [17]  on  this  graph.  That 
is,  each  vertex  picks  a  random  number  between  1  and  j ,  where  j  is  large  enough  that  the 
probability  of  ties  is  small.  If  v  has  a  number  higher  than  any  of  its  neighbors’,  it  enters  the 
MIS  and  its  neighbors  drop  out.  We  then  iterate. 

Let  M  C  V'  denote  the  resulting  Mis.  For  any  v  G  V7,  C*  intersects  C*  internally  for  at 
most  24  other  vertices  u  G  V'\  thus  with  probability  A  _  o(l),  v  will  enter  M  on  the  first 
iteration.  Moreover,  if  u,v  G  V'  are  at  a  distance  of  at  least  1 17  log  n  from  each  other,  then 
they  have  no  common  neighbors  in  G',  and  so  these  events  are  independent.  Thus, 

Lemma  3.2  Let  u,v  G  V'  be  such  that  d(u,v )  >  llylogn.  Then  with  constant  probability, 
both  u  and  v  belong  to  the  set  M  constructed  above. 

If  v  G  M,  we  will  call  Cv  a  cluster.  We  now  want  to  construct  internally  disjoint  enclosures 
Dv  around  each  C*,  for  v  G  M,  such  that  every  vertex  of  G  belongs  to  some  enclosure  and 
such  that  each  Dv  is  a  union  of  7-blocks.  The  sets  C*  are  disjoint  and  are  unions  of  7-blocks, 
but  they  do  not  cover  all  of  G.  However,  by  the  maximality  of  M,  any  7-block  X  that  does 
not  belong  to  C*  for  some  v  G  M  must  share  a  boundary  vertex  with  such  a  C*.  For  each 
such  X,  we  pick  such  a  C*  arbitrarily  and  add  X  to  Dv.  Thus  the  Dv  now  form  a  partition 
of  G,  and  each  Dv  is  a  union  of  7-blocks. 

We  now  define  J\f  to  be  the  graph  on  vertex  set  M,  with  u  and  v  joined  by  an  edge  if 
some  7-block  of  Du  shares  a  wall  with  a  7-block  of  Dv.  As  argued  above,  any  7-block  in  Dv 
must  belong  to  the  5  x  5  set  of  7-blocks  centered  at  Gm  and  from  this  it  is  easy  to  argue 
that  at  most  20  7-blocks  not  contained  in  Dv  can  share  a  wall  with  a  7-block  of  Dv.  Thus 
we  have 

Lemma  3.3  The  degree  of  a  vertex  in  A f  is  at  most  20. 

If  we  were  to  contract  each  enclosure  Dv ,  the  resulting  graph  would  contain  all  the  edges 
of  J\f  with  multiplicity  0 ( log  n)  —  enough  to  run  a  “high-bandwidth”  algorithm.  But  given 
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a  routing  in  A/",  we  then  run  into  the  following  problem:  we  need  to  convert  the  sequence  of 
neighboring  clusters  connecting  the  clusters  of  st  and  t,  in  the  contracted  graph  back  to  a 
routing  in  G.  For  this  we  use  the  natural  “crossbar”  structures  in  the  mesh  in  the  on-line 
case,  with  additional  flow  techniques  in  the  off-line  case. 

We  start  by  developing  the  “crossbar”  structures  we  use. 

Definition  3.4  A  v-ring  is  a  subgraph  G[X\,  where  X  is  the  set  of  vertices  at  distance 
exactly  r  from  v,  for  some  r  between  7^7  log  n  and  ylogn.  Thus  a  v-ring  is  either  a  cycle  or 
a  path,  depending  on  whether  Cv  has  a  wall  on  the  boundary  of  G.  If  R  and  R'  are  v-rings, 
we  say  that  R  is  inside  (resp.  outside,)  R'  if  the  distance  from,  R  to  v  is  less  than  (resp. 
greater  than)  the  distance  from  R'  to  v. 

Note  that  the  set  of  v-rings  are  the  log  n  disjoint  cycles  around  v  right  outside  the  bound¬ 
ary  of  Cv,  in  the  “inner  half”  of  C*  \CV. 

At  this  point  we  will  need  some  additional  notation:  If  A"  C  V,  let  <5(X)  denote  the  set 
of  edges  leaving  S,  and  vr(A")  the  set  of  vertices  of  X  incident  to  an  edge  of  h(A").  For  two 
sets  X,  Y  C  V,  let  <5(X,  Y)  denote  the  set  of  edges  crossing  from  X  to  Y . 

For  each  pair  ( v ,  w )  that  is  an  edge  of  A I,  we  choose,  for  a  sufficiently  small  constant  p, 
a  set  tV:W  of  plogn  edges  in  5(DV,  Dw ).  Let  t'v  denote  the  set  of  all  vertices  in  Dv  that  are 
incident  to  an  edge  in  some  tv,w.  We  also  choose  a  set  a'v  of  plogn  vertices  evenly  spaced 
on  the  outer  boundary  of  Cv.  Now  by  Lemma  3.3,  we  can  choose  p  small  enough  that 
| t'v  U  cr'|  <  lylogn,  and  hence  we  can  associate  a  different  v-ring  to  each  vertex  in  t'v  U  cr'v. 
Moreover,  in  a  straightforward  fashion  we  can  construct  edge-disjoint  paths  from  each  such 
vertex  to  its  associated  ring.  We  assign  the  outermost  plog  n  rings  to  the  vertices  of  a'v.  For 
u  G  r'  U  a'v,  let  Y)f  denote  the  union  of  the  ring  associated  with  u  with  the  path  from  u  to 
this  ring.  Then  we  have 

Lemma  3.5  The  (non-simple)  paths  Yjf  are  mutually  edge- disjoint,  and  each  pair  meets  at 
some  vertex  of  C*  \  Cv. 

Proof.  The  paths  are  edge-disjoint  by  construction.  Now  suppose  u,w  G  r,',  and  that  the 
ring  associated  with  u  is  inside  the  ring  associated  with  w.  Then  the  path  from  u  to  its  ring 
must  intersect  the  ring  associated  with  w,  and  so  Y)f  and  Y™  meet  at  a  vertex.  The  same 
argument  applies  if  u,w  G  <r',,  and  also  if  u  G  cf'v  and  tv  G  t'v  because  in  this  case  the  ring 
associated  with  u  lies  outside  the  ring  associated  with  tv.  ■ 

We  are  now  ready  to  describe  and  analyze  the  routing  algorithms  themselves. 

3.2  The  On-Line  Algorithm 

Say  that  a  terminal  pair  ( Si,U )  G  T  is  short  if  dfsi,  A)  <  lGylogn,  and  long  otherwise.  The 
on-line  algorithm  makes  an  initial  random  decision  whether  to  accept  only  short  connections 
or  only  long  connections;  this  costs  at  most  a  factor  of  two  in  the  competitive  ratio.  Below 
we  give  0(log  n)-competitive  algorithms  for  handling  each  type  of  connection. 
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3.2.1  Routing  Long  Connections 

First,  we  only  consider  terminal  pairs  with  both  ends  in  sets  of  the  form  Cv  —  denote  this 
set  of  terminal  pairs  by  Tm ■  If  T*  denotes  a  realizable  subset  of  maximum  size,  then  by 
Lemma  3.2,  the  expected  number  of  terminal  pairs  in  T*  that  belong  to  7m  is  a  constant 
fraction  of  \T*\.  Thus  we  only  lose  a  constant  factor  in  the  competitive  ratio  by  restricting 
attention  to  TM. 

Let  A f{c)  denote  the  graph  J\f  in  which  each  edge  has  been  given  a  capacity  of  c.  We 
now  define  an  on-line  routing  problem  in  the  simulated  network  J\f (p  log  n).  If  s*  G  Cv, 
then  we  define  its  image  in  the  “simulation”  to  be  'ijj(si)  =  v.  The  input  will  simply  be  the 
sequence  of  terminal  pairs  (/ip(si),'ip(ti)),  where  ( Si,ti )  is  the  sequence  of  pairs  presented  to 
the  algorithm  running  on  G.  Our  algorithm  for  the  problem  in  the  simulated  network  is  as 
follow:  we  route  (v,w)  if  (i)  the  AAP  algorithm  on  Af  (plogn)  accepts  (v,w),  and  (ii)  no 
connection  with  an  end  equal  to  either  v  or  w  has  yet  been  accepted. 

Lemma  3.6  The  above  algorithm  is  O (log n)- competitive  against  the  fractional  optimum  in 
M(p  log  n) . 

Proof.  Let  (v,  w)  be  a  requested  connection  in  J\f .  Following  the  approach  developed  in 
Section  2.2,  we  view  ( v ,  w)  as  being  “judged”  by  three  on-line  algorithms:  the  AAP  algorithm, 
an  algorithm  Av  which  only  permits  one  connection  with  an  end  equal  to  v,  and  an  algorithm 
Aw  which  only  permits  one  connection  with  an  end  equal  to  w.  By  Proposition  2.1,  the  AAP 
algorithm  is  0(logn)-competitive;  the  algorithms  Av  and  Aw  are  also  O ( log  n)-eornpetitive 
since  the  maximum  fractional  weight  of  connections  that  the  optimum  can  accept  originating 
at  any  one  node  of  J\f  is  O(logn).  Thus,  applying  Proposition  2.3  with  c*  =  O(logn)  and 
d  —  3,  we  see  that  the  combined  routing  algorithm  is  0(logn)-competitive  against  the 
fractional  optimum.  ■ 

Our  on-line  algorithm  in  G  simply  runs  the  above  simulation;  whenever  ('if(si),  'f’itf))  is 
accepted,  it  routes  the  pair  (s*,  U)  in  G  using  the  paths  constructed  in  Lemma  3.5.  The 
following  lemma  says  that  it  will  not  run  out  of  “bandwidth”  while  doing  this. 

Lemma  3.7  The  algorithm  in  G  can  route  all  the  connections  accepted  by  the  simulation. 

Proof.  When  the  simulation  accepts  (ip(si),  "0(A)),  it  specifies  a  sequence  of  neighboring 
clusters  CVl,  CV2, . . . ,  CVr,  where  v\  =  'if(si)  and  vr  =  if{U). 

The  algorithm  in  G  routes  by  simply  moving  from  one  cluster  in  this  sequence 

to  the  next  using  paths  of  the  form  Yf;.  More  concretely,  it  first  chooses  any  w  G  a'vi  and 
w'  G  a’  and  sets  Z0  =  Y™  and  Zr  =  Y™\  Then  for  each  j  =  1, ...  ,  r  —  1,  it  chooses  any 
edge  (u,u')  G  rVjtVj+1  that  has  not  yet  been  used  for  routing,  and  sets  Zj  =  Yff  UY™'+1.  Since 
there  are  at  least  plogn  such  edges  available,  and  the  simulation  accepts  at  most  plogn  pairs 
whose  routes  use  the  edge  in  J\f  from  Vj  to  Vj+ 1,  this  is  always  possible.  Now,  by  Lemma  3.5, 
the  union  Z0  U  •  •  •  U  Zr  contains  a  path  from  the  boundary  of  CVl  to  the  boundary  of  CVr.  ■ 

Finally,  we  have  to  show  that  optimum  in  G  is  not  far  from  the  optimum  in  the  simulation. 
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Lemma  3.8  For  any  realizable  subset  T  ofTM,  if( T ')  can  be  routed  in  J\f(5^/\ogn). 

Proof.  For  each  path  P  in  the  optimal  routing,  construct  the  following  path  for 
^(ti))  in  J\f  —  when  P  crosses  from  Dw  into  Dw>,  add  an  edge  from  w  to  w' .  Since 
\S(DW,  Dwf)\  <  57  log n,  so  at  most  this  many  paths  in  the  constructed  routing  will  use  the 
edge  (w,  w ')  in  AT.  m 

Now,  since  the  on-line  algorithm  is  0(logn)-competitive  against  the  fractional  optimum 
in  Af(plogn),  it  is  also  0(logn)-competitive  against  the  fractional  optimum  in  A/”(57logn), 
which  by  Lemmas  3.2  and  3.8  is  at  least  as  large  as  the  maximum  realizable  subset  of  T. 
Thus  we  are  0(logn)-competitive  in  routing  long  connections. 

з. 2.2  Routing  Short  Connections 

To  handle  short  connections,  we  require  the  following  two  facts. 

Proposition  3.9  Let  H  =  (V,  E )  be  an  arbitrary  graph  of  diameter  d.  Then  there  is  a 
deterministic  on-line  mdpp  algorithm  that  is  2  ■  max(d,  \J\E\ )  -  competitive. 

Proof.  Let  m  =  \E\.  The  algorithm  maintains  a  sequence  of  graphs  H\,  H-2, . . .  as  follows. 
H\  =  H.  The  algorithm  always  routes  the  first  request  on  a  shortest  path  Pi,  and  sets 
H-2  —  H\  \P\.  In  general,  when  presented  with  request  ( Si,U ),  the  algorithm  routes  it  on  a 
shortest  path  Pj  in  Hi  if  dfsi,ti)  <  \Jm  in  H;.  It  then  sets  Hi+l  =  Hi  \  Pi.  Let  p  denote  the 
total  number  of  paths  routed  by  the  algorithm. 

Let  d'  =  max(d,  y/rri).  Consider  any  routing  for  T,  consisting  of  paths  Qi, ... ,  Qq.  At 
most  d'  of  the  Q}  intersect  each  Pi,  since  the  Qi  are  all  edge-disjoint  and  jP,;|  <  d' .  Also, 
at  most  yjm  of  the  Qj  fail  to  intersect  any  of  the  P,;,  since  the  pair  ( Sj,tj )  associated  with 
Qj  must  have  been  rejected  by  the  on-line  algorithm,  and  hence  \Qj\  >  \fm.  Thus  we  have 
q  <  d'p  +  yfm  <  2 d'p.  ■ 

Lemma  3.10  Let  X  =  G[{i  —  r )  :  ( i  +  r),  (j  —  r)  :  (j  +  r)]  for  some  natural  numbers  i,  j, 
and  r;  let  Y  =  G[{i  —  2 r)  :  ( i  +  2r),  (j  —  2 r)  :  (j  +  2r)];  and  let  T  be  a  set  of  terminal  pairs 
with  both  ends  in  X .  Then  the  maximum  size  of  a  subset  of  T'  that  is  realizable  in  G\Y]  is 
at  least  |  the  maximum  size  of  a  subset  of  T  that  is  realizable  in  G. 

Proof.  Let  T"  C  T  denote  a  realizable  set  of  maximum  size,  and  consider  the  set  of  paths 
in  a  routing  of  T" .  Consider  the  set  of  paths  leaving  X.  Since  |h(Af)|  <  8rj  there  are  at 
most  4 r  such  paths.  Delete  the  portion  of  each  path  between  its  first  and  last  intersection 
with  h(Af);  using  the  r  disjoint  rings  in  G(V\  AT]  (as  defined  in  the  previous  section),  we  can 
connect  the  resulting  “breakpoints”  of  at  least  r  of  these  paths  using  edge-disjoint  paths  in 
G[Y  \X],  Thus,  we  are  still  routing  at  least  |  of  T"  in  G\Y].  ■ 

The  algorithm  for  short  connections  is  now  as  follows.  Set  r  =  32y  log  n,  and  run  Luby’s 
algorithm  to  find  a  subset  M'  of  V,  maximal  subject  to  the  property  that  Loo(u,u)  >  4r  for 

и, v  G  M' .  For  u  G  M',  dehne  Xu  to  be  the  set  of  all  vertices  whose  L0 0  distance  from  u  is 
at  most  r,  and  Yu  to  be  the  set  of  all  vertices  whose  L ^  distance  from  u  is  at  most  2 r. 
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Lemma  3.11  If  (si,ti)  is  a  short  connection,  then  with  constant  probability  there  is  an  such 
that  Si,ti  G  Xu. 

Proof.  A  sufficient  condition  for  some  Xu  containing  both  sl  and  ti  to  enter  the  MIS  in  the 
first  iteration  is  that  among  all  vertices  within  L0 0  distance  4r  of  Sj,  the  one  that  picks  the 
highest  number  has  distance  at  most  |r  from  st.  This  happens  with  constant  probability. 


Let  %,  denote  the  set  of  short  connections  both  of  whose  ends  lie  in  Xu.  We  now  run 
the  algorithm  of  Proposition  3.9  on  the  (disjoint)  subgraphs  G \YV]  simultaneously,  using 
the  Tu  as  the  sets  of  terminal  pairs.  By  Proposition  3.9  and  Lemma  3.10,  we  are  O(logn)- 
competitive  in  each  subgraph.  Thus,  by  Lemma  3.11,  we  are  O (log  n)-competitive  in  routing 
short  connections.  Thus, 

Theorem  3.12  The  on-line  algorithm  is  O (log  n)- competitive  in  the  two-dimensional  mesh. 


3.3  The  Off-Line  Algorithm 


For  the  constant-factor  off-line  approximation,  we  use  a  variant  of  the  graph  A/".  In  J\f( 7  log  n) , 
for  any  fixed  constant  7,  one  can  obtain  a  constant-factor  approximation  to  the  MDPP  by  the 
following  randomized  rounding  algorithm  of  Raghavan  and  Thompson  [23,  22],  First  we  solve 
the  fractional  relaxation  of  the  MDPP  instance  (this  can  be  done  in  polynomial  time);  from 
this,  we  obtain  for  each  terminal  pair  (sj,  ti)  a  collection  of  paths  Pf, . . . ,  Pf  and  associated 
weights  y\, . . . ,  yf  G  [0, 1]  such  that  Xi  =  Jfj  y\  G  [0, 1].  We  now  pick  a  scaling  factor  y  <  1; 
independently  for  each  terminal  pair  (s*,  tf)  we  route  it  on  path  P-  with  probability  yyj,  and 
don’t  route  it  at  all  with  probability  1  —  yxi.  If  we  do  route  it,  we  say  that  (sj,fj)  has  been 
rounded  up.  In  [23,  22],  it  is  shown  that  with  constant  probability,  no  capacity  is  violated 
by  the  selected  paths,  and  the  number  of  pairs  that  are  rounded  up  is  a  constant  fraction  of 
the  fractional  optimum. 

In  particular,  we  require  the  following  theorem  from  [22], 


Theorem  3.13  (Raghavan)  Let  X \ ,  X2 ,  ■  ■  • ,  Xr  be  independent  Bernoulli  trials  with  EXj  = 
Pj  and  T  =  X so  E^>  =  m  —  ff  iPi-  Then  for  5  >  0  we  have 


Pr  [T  >  (1  +  S)m  < 


es 

(1  +  <5)d+5) 


We  specialize  this  to  the  form  in  which  we  will  use  it  as  follows. 

Corollary  3.14  Let  0  <  y  <  1,  and  pi,  ■  ■  ■  ,pr  G  [0, 1].  Let  X[,X'2, . . . , X'r  be  independent 
Bernoulli  trials  with  EX^  =  ypj  and  T'  =  t  X) .  Let  m  =  Jf  iPi-  Then 

Pr[T;  >  m]  <  (ey)m. 
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Proof.  Apply  the  bound  of  Theorem  3.13  with  m  set  to  pm  and  <5  set  to  /i  1  —  1.  ■ 

Let  us  consider  how  to  use  this  randomized  rounding  approach  in  routing  long  connections 
off-line.  In  the  high-capacity  network  A/",  this  rounding  approach  is  fine;  but  to  get  a  constant- 
factor  approximation  we  also  have  to  be  within  a  constant  factor  of  the  optimum  in  routing 
terminals  out  of  the  clusters  (in  the  on-line  algorithm  it  was  enough  to  route  only  one). 
To  this  end,  we  build  the  following  more  complicated  network  J\P.  Let  zv  denote  the  node 
representing  v  G  M  in  the  network  A f\  we  construct  A/7  by  attaching  Cv  to  zv  via  an  edge 
from  zv  to  each  node  in  n(Cv).  Let  A/7 (7)  denote  the  network  A/7  in  which  each  edge  between 
nodes  of  the  subgraph  J\f  has  capacity  7,  and  all  other  edges  have  unit  capacity. 

Let  f*  denote  the  value  of  the  fractional  optimum  solution  to  the  mdpp  in  Af'(plogn). 
We  now  run  the  randomized  rounding  algorithm  on  Af(plogn).  With  high  probability  the 
parts  of  all  selected  paths  lying  in  the  subgraph  Af(plogn),  taken  together,  do  not  violate 
any  capacity  constraint;  and  the  number  of  pairs  that  are  rounded  up  is  within  a  constant 
factor  of  the  fractional  optimum  f*.  We  now  must  convert  the  selected  paths  in  Af'(plogn) 
into  Si-ti  paths  in  G.  We  can  use  the  technique  of  the  previous  section  to  produce,  for 
each  selected  pair  (s^tf),  a  “global”  path  Pj  that  begins  at  rG  ^  C  and  ends  at 

TV>(t;)  G.  TT 

The  real  problem  is  how  to  find  paths  within  the  clusters  such  that  each  Si  (resp.  A) 
that  has  been  rounded  up  can  reach  the  endpoint  of  this  associated  global  path  on  tLs.^ 
(resp.  rW.)).  For  this,  the  paths  returned  by  the  randomized  rounding  are  of  no  value,  since 
the  edges  of  A/7  within  the  clusters  Cv  have  only  unit  capacity.  Instead  we  argue  as  follows. 

Let  Sv  denote  the  set  of  terminals  in  Cv  that  are  rounded  up.  Each  is  trying  to  “get 
out”  to  its  associated  path  that  begins  at  r'.  Given  the  crossbar  in  C*  \CV,  it  is  sufficient 
to  route  each  terminal  in  Sv  to  any  vertex  in  a'v. 

So  this  leaves  us  with  the  following  escape  problem.  We  are  given  the  set  Sv  of  terminals 
that  have  been  rounded  up,  and  we  want  to  route  a  large  fraction  of  them  to  a'v.  The 
following  lemma,  whose  proof  contains  the  central  step  of  the  algorithm,  says  that  this  can 
be  done. 

Lemma  3.15  For  a  sufficiently  small  (constant)  value  of  p,  there  is  a  constant  c  <  1  so 
that  we  can  find  sets  S'v  C  Sv,  such  that 

(i)  if  one  end  of  a  pair  (si,U)  belongs  to  then  so  does  the  other, 

(ii)  each  set  S'v  can  be  linked  to  a'v  C  tt (Cv)  via  edge-disjoint  paths. 

(Hi)  the  probability  that  U,,5'(|  >  cf*  is  at  least  a  constant,  where  the  probability  is  taken 
over  the  randomized  rounding  that  produced  the  sets  Sv. 

Proof.  We  will  first  construct  such  a  set  with  condition  (iii)  weakened  to  the  requirement 
that  each  S'v  can  be  linked  to  n(Cv)  via  edge-disjoint  paths  (rather  than  to  the  smaller  set 
a'y).  This  is  sufficient  to  imply  the  lemma,  as  follows.  Suppose  we  obtain  such  sets  S".  For 
each  s  G  S ",  identify  the  vertex  in  a'v  closest  to  s.  We  now  build  a  graph  on  the  set  of 
terminal  pairs  in  U VS",  joining  two  if  at  either  end  they  have  the  same  closest  vertex  in  some 
set  a'v.  We  claim  this  graph  has  degree  at  most  87p_1;  for  the  spacing  between  vertices  of 
a'v  on  the  boundary  of  Cv  is  47/?-1,  and  so  at  most  twice  this  number  of  terminal  pairs  can 
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compete  with  some  (at  either  end)  for  the  same  vertex  of  a'v.  Thus  this  graph  has 

an  independent  set  of  size  at  least  8  i1+1  |Ut>)S,^|;  if  we  let  the  terminal  pairs  in  this  large 
independent  set  constitute  the  sets  Sv,  we  satisfy  the  requirements  of  the  lemma. 

This  allows  us  to  deal  with  a  standard  escape  problem  on  a  rectangular  mesh:  each 
terminal  is  allowed  to  “escape”  to  any  vertex  on  the  boundary.  First  observe  the  following 
fact:  an  escape  problem  on  a  rectangular  mesh  is  feasible  if  and  only  if.  for  all  p,  q,  any 
subrectangle  of  size  p  x  q  contains  at  most  2 (p  +  q)  terminals.  To  see  this,  note  that  we 
can  reduce  the  escape  problem  to  a  single-source/single-sink  maximum  flow  problem,  and 
thus  only  have  to  verify  the  cut  condition.  On  a  rectangular  mesh,  the  smallest  rectangle 
enclosing  any  connected  cut  has  no  greater  capacity,  and  contains  at  least  as  many  terminals, 
as  the  original  cut;  thus  the  cut  condition  holds  if  and  only  if  it  holds  for  all  subrectangles. 

This  suggests  the  following  algorithm  to  construct  the  set  S".  Call  a  rectangle  overfull 
if  it  violates  the  cut  condition;  we  go  through  each  s  e  Sv ,  deleting  it  if  it  is  contained  in 
any  overfull  rectangle  —  we  also  then  delete  its  matching  terminal  in  some  other  cluster. 
This  results  in  a  set  S",  which  by  the  argument  of  the  previous  paragraph  can  be  completely 
routed  to  tt(Cv)  on  edge-disjoint  paths. 

ft  remains  to  lower-bound  the  expected  size  of  LbS1".  Say  that  a  terminal  s  survives  if 

(i)  it  is  initially  rounded  up,  and  hence  included  in  a  set  Sv,  and 

(ii)  it  is  not  deleted  because  it  or  its  matching  terminal  is  contained  in  an  overfull  rect¬ 
angle. 

So  what  is  the  probability  that  a  terminal  s,  which  is  given  weight  fs  G  [0, 1]  by  the 
fractional  optimum,  survives?  First  we  consider  the  probability  that  s  is  contained  in  a  fixed 
overfull  px  q  rectangle  R ,  given  that  s  has  been  rounded  up.  In  order  for  this  rectangle  R  to 
become  overfull,  it  must  be  that  at  least  2 (p  +  q)  terminals  in  R  other  than  s  were  rounded 
up.  But  since  the  un-scaled  fractional  flow  was  feasible,  the  total  fractional  weight  contained 
in  R  is  at  most  2p(p  +  q).  Thus,  setting  y  =  (ep)2,  the  probability  that  the  rectangle  R 
becomes  overfull  after  rounding,  given  that  s  is  rounded  up,  is  at  most  yp+q. 

Now,  since  s  is  contained  in  at  most  pq  rectangles  of  dimensions  p  x  q,  the  probability 
that  s  is  contained  in  any  overfull  rectangle  after  rounding,  given  that  it  is  rounded  up,  is 
clearly  bounded  by  the  infinite  sum 


J2J2,p(iyp+q  = 

■p  q 


(1-2/)4' 


s  can  also  be  deleted  if  its  matching  terminal  is  contained  in  an  overfull  rectangle,  so  the 
probability  of  s  being  deleted,  after  having  been  rounded  up,  is  at  most  jpzyp  •  By  taking  p 
small  enough,  we  can  make  this  last  expression  a  constant  less  than  . 

Finally,  the  probability  that  s  survives  is  equal  to  the  probability  that  it  is  rounded  up, 
which  is  pfs,  times  the  probability  that  it  is  not  deleted  after  being  rounded  up,  which  is  at 
least  Thus,  the  expected  size  of  UvS'f  is  at  least  \ pf  *,  for  a  small  enough  constant  value 
of  p.  So  by  Markov’s  inequality,  the  probability  that  |U1,,S''/|  >  \pf*  is  at  least  an  absolute 
constant;  and  the  lemma  follows.  ■ 
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To  turn  the  above  lemma  from  a  statement  holding  with  constant  probability  to  one 
with  high  probability,  we  repeat  the  randomized  rounding  0(\og  n)  times  and  take  the  best 
routing  obtained.  A  single  run  of  randomized  rounding  fails  if  one  of  the  following  two  bad 
events  happens:  (1)  one  of  the  high  capacity  edges  of  the  simulated  network  is  used  by  too 
many  selected  paths,  or  (2)  the  set  LbA'  selected  by  Lemma  3.15  is  too  small.  The  first  event 
has  inverse  polynomially  small  probability,  and  the  probability  of  the  second  is  bounded  by 
a  constant,  so  the  total  probability  of  any  bad  event  is  bounded  by  a  constant  less  than  1. 
This  implies  that  out  of  the  O(logn)  tries,  one  must  succeed  with  high  probability. 

This  gives  a  constant-factor  approximation  for  long  connections:  if  (s;,fj)  is  rounded  up, 
and  Si,ti  G  U„S'',  then  we  concatenate  the  paths  from  Sj  to  7r(C^(Si))  (given  by  Lemma  3.15) 
to  7r(D^(Sj))  (given  by  the  crossbar  in  C*  \  Cv)  to  n{D^tif)  (given  by  the  edges  of  the  path 
in  A/"(plogn),  joined  together  with  crossbars  as  in  Lemma  3.7),  and  now  symmetrically  to 
n(C^ti))  and  to  f. 

We  handle  short  connections  as  in  Section  3.2.2.  That  is,  we  use  a  randomized  algorithm 
to  construct  subsquares  Xu  C  Yu,  with  all  Yu  disjoint,  and  only  handle  short  connections 
both  of  whose  ends  lie  in  a  single  Xu.  In  this  case,  however,  we  now  run  the  above  algorithm 
recursively  on  each  G\Yy\. 

Call  a  connection  “medium”  if  it  is  now  a  long  connection  in  this  recursive  call,  and 
“small”  otherwise.  Medium  connections  are  handled  as  described  above.  Small  connections 
take  place  within  clusters  of  size  O (log  log  n)  and  therefore  can  be  simply  solved  to  optimality 
by  brute  force.  We  can  then  take  the  largest  realizable  set  we  find  among  the  long,  medium, 
and  small  connections,  obtaining 

Theorem  3.16  There  is  a  randomized  (off-line)  mdpp  algorithm  in  the  two-dimensional 
mesh  that  produces  a  constant-factor  approximation  with  high  probability. 

4  Densely  Embedded  Graphs:  Definition  and  Proper¬ 
ties 

We  now  want  to  extend  the  above  algorithm  to  any  graph  that  is  sufficiently  “mesh-like.”  We 
define  the  class  of  graphs  here;  in  the  following  section  we  show  how  to  extend  the  routing 
algorithms  to  this  class. 

We  will  need  some  additional  notation  and  definitions:  Recall  that  for  u,  v  G  V,  d(u,  v) 
denotes  the  least  number  of  edges  in  a  u-v  path.  Let  Br{v)  =  {u  :  d(u,v )  <  r}.  Also  recall 
that  <5(A")  and  vr(A")  denote  the  set  of  edges  leaving  A"  and  the  set  of  vertices  of  X  incident  to 
edges  in  <5(A").  We  write  X°  —  S  \  vr(A").  Observe  that  removing  vr(X)  from  X  disconnects 
it  from  the  rest  of  the  graph,  and  n(Br(v))  consists  of  vertices  at  exactly  distance  r  from 
v.  If  C  is  a  connected  subset  of  G  \  X,  we  use  T(X,  C)  to  denote  the  (unique)  connected 
component  of  G  \  X  containing  C.  The  set  of  vertices  in  vr(X)  which  have  a  neighbor  in 
T(X,  C)  will  be  called  the  segment  of  tt(X)  bordering  C  and  denoted  a(X,  C).  We  say  that 
a  set  X  C  V  is  simple  if  G  \  X  is  connected. 
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Definition  4.1  A  graph  H  is  an  a- semi- expander  if  for  every  X  C  V(H)  for  which  |A"j  < 
\\V(H)\,  we  have  |<5(A)|  >  ot\J\X\. 

Since  our  goal  is  to  generalize  the  two-dimensional  mesh,  let  us  note  the  following  prop¬ 
erties  of  the  mesh. 

(i)  It  is  a  planar  graph  with  bounded  degree,  and  (aside  from  one  “exceptional  face” )  it 
is  Eulerian  and  has  bounded  face  size. 

(ii)  It  is  an  o-semi-expander,  for  a  constant  a  >  0  based  on  the  ratio  of  the  two  side 
lengths  of  the  mesh. 

(iii)  Square  sub- meshes  of  the  mesh  satisfy  (i)  and  (ii). 

In  the  arguments  to  follow,  it  is  quite  cumbersome  —  though  not  technically  difficult  - 
to  deal  with  “exceptional  faces”  of  the  type  in  (i).  Thus,  in  the  following  section  we  will 
work  with  the  more  restricted  class  of  uniformly  densely  embedded  graphs ,  where  all  faces 
have  bounded  size;  and  we  will  assume  further  assume  that  G  is  Eulerian.  In  Section  6.2, 
we  show  how  to  handle  graphs  with  an  exceptional  face;  in  this  way,  our  class  of  graphs  will 
include  the  two-dimensional  mesh. 

First  we  need  some  preliminary  topological  definitions.  Let  E  denote  a  compact  orientable 
surface;  it  is  well-known  (see  e.g.  [18])  that  E  may  be  obtained  from  the  2-sphere  by  attaching 
a  finite  number  of  handles.  By  a  disc  we  mean  a  set  homeomorphic  to  the  closed  unit  ball  in 
R2,  and  by  E -disc,  we  mean  a  subset  of  E  homeomorphic  to  a  disc.  Our  definition  of  graph 
embedding  is  standard;  a  face  of  an  embedded  graph  G  is  a  connected  component  of  E  \  G, 
and  we  say  G  is  strongly  embedded  if  the  closure  of  each  face  is  a  E-disc,  and  each  face  is 
bounded  by  a  simple  cycle  of  G.  Finally,  we  also  use  the  the  terms  curve  (continuous  image 
of  [0, 1])  and  closed  curve  (continuous  image  of  S1). 

Our  class  of  graphs  is  defined  to  satisfy  analogues  of  properties  (i),  (ii),  and  (iii)  locally. 

Definition  4.2  A  graph  G  =  (V,  E)  is  uniformly  densely  embedded  with  parameters  a,  A, 
A,  and  £  if: 

(i)  G  is  strongly  embedded  on  a  compact  orientable  surface  E,  it  has  maximum  degree  A, 
and  each  face  is  bounded  by  at  most  £  edges. 

(ii)  For  each  r  <  A  log n  and  each  v  G  V,  the  drawing  of  G[Br(v)]  is  contained  in  a 
'E-disc. 

(iii)  For  each  r  <  A  log  n  and  each  v  G  V,  the  graph  G[Br(v)\  is  an  a- semi- expander. 

Thus,  for  the  remainder  of  the  paper  aside  from  Section  6.2,  we  will  assume  that  G  is  a 
simple  Eulerian  graph  that  is  uniformly  densely  embedded  on  a  surface  E  with  parameters  a, 
A,  A,  and  t.  In  Section  6.2,  we  show  how  our  algorithms  can  be  adapted  to  handle  graphs 
satisfying  the  following  weaker  definition;  it  is  the  same  as  the  definition  above,  except  that 
we  allow  an  exceptional  face. 

Definition  4.3  A  graph  G  =  (V,  E )  is  densely  embedded  and  nearly-Eulerian  with  param¬ 
eters  a,  A,  A,  and  £  if: 

(i)  G  is  strongly  embedded  on  a  compact  orientable  surface  E  and  has  maximum  degree 

A. 
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(%)'  G  contains  a  face  $*  such  that  all  faces  other  than  are  bounded  by  at  most  £  edges, 
and  every  vertex  ?iot  incident  to  <F*  has  even  degree. 

(ii)  For  each  r  <  Alogn  and  each  v  G  V,  the  drawing  of  G[Br(v)]  is  contained  in  a 
E -disc. 

(in)  For  each  r  <  Alogn  and  each  v  G  V ,  the  graph  G[Br{v)\  is  an  a-semi-expander. 

The  classes  of  graphs  satisfying  these  definitions  are  incomparable  to  the  class  considered 
in  our  earlier  paper  [14].  The  semi- expansion  condition  above  will  be  shown  to  imply  the 
uniformly  high- diameter  condition  of  [14]  (see  Lemma  4.4);  however,  in  the  current  paper,  we 
only  require  planarity  and  semi-expansion  locally,  and  essentially  no  restrictions  are  placed 
here  on  the  “global”  structure  of  the  graph.  The  examples  of  uniformly  high-diameter  graphs 
constructed  in  [14]  are  densely  embedded  as  well;  and  in  Section  4.2  we  will  discuss  some 
related  classes  of  graphs  that  are  densely  embedded. 

4.1  Some  Basic  Properties 

We  now  show  that  our  definition  implies  G  has  some  additional  mesh-like  properties.  First  of 
all,  for  any  v  G  V  and  r  <  Alogn,  the  fact  that  G[Br(v)\  is  a  bounded-degree  semi-expander 
implies  that  the  set  Br{y )  has  size  at  least  quadratic  in  r;  by  also  using  the  planarity  of 
G[Br(v)\,  one  can  show  an  analogous  upper  bound.  We  summarize  this  as  follows. 

Lemma  4.4  There  are  constants  a  and  /3  depending  only  on  a  and  A  such  that  the  following 
holds.  For  each  r  <  Alogn  and  each  v  G  V,  we  have  ar 2  <  \Br(v)\  <  (3r 2 . 

Proof.  Fix  r  <  Alogn  and  v  G  V,  and  let  S  =  Br{y).  To  see  the  lower  bound,  note  that  for 
any  %  <  r,  if  Xi  —  |B,:(u)|,  then  by  the  semi-expansion  of  F[  we  have 

^  a  , _ 

Xi  >  Xi-!  +  —  _  y/Xi- 1- 

For  at  least  a^/Xi-i  edges  leave  Bi-i(v),  and  at  most  A  —  1  are  incident  to  any  one  vertex. 
Let  v  =  \  now  one  verifies  by  induction  that  x%  > 

Xi  >  -  l)2  +  (n2(i  -  1) 

=  L^(i-l)(i  +  3) 

1 

-  L  2  -2 
>  - V  l  . 

~  16 

To  see  the  upper  bound,  we  observe  that  GtS1]  is  planar  and  has  diameter  at  most  2 r. 
Let  n  =  \S\.  By  the  Lipton-Tarjan  planar  separator  theorem  [16],  there  is  a  set  of  at  most 
4r  +  1  vertices  whose  removal  breaks  H  into  components  each  of  size  at  most  | n.  Let  X  be 
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a  union  of  these  components  of  size  between  and  |n.  Then 


<  IWI 

<  (A  —  l)(4r  +  1) 

from  which  the  result  follows.  ■ 

The  following  two  facts  are  quite  useful;  the  first  essentially  relates  the  size  of  the  “perime¬ 
ter”  of  a  set  Br(y)  ( r  <  A  logn)  to  its  radius. 

Lemma  4.5  Let  c  >  1  and  r  a  positive  integer  be  such  that  cr  <  A  log  n.  Then  for  some  r' 
between  r  and  cr,  we  have  \ir(Bri(v))\  <  (3  ■  ■  r. 

Proof.  Since  n (Br>(v))  C  {u  :  d(v,u )  =  r'},  the  sets  n(Br(v)),  n(Br+i(v)), . . . ,  /n(Bcr{v)) 
are  all  disjoint  and  contained  in  Bcr{y).  Since  \Bcr(y)\  <  (3c2,r2,  one  of  these  sets  has  size  at 
most  (3  ■  ■  r.  m 

Using  this,  we  show  that  we  can  extend  any  small  enough  set  U  to  a  simple  set  with  at 
most  a  constant-factor  increase  in  its  radius.  Recall  from  the  beginning  of  Section  4  that  a 
set  X  is  simple  if  G  \  X  is  connected. 

Lemma  4.6  There  is  a  constant  f  such  that  the  following  holds.  Let  U  C  Br(y),  where 
r  <  ^A  logn.  Then  there  is  a  component  T  of  G\U  and  a  planar  simple  set  U'  such  that 
f/cl/'c  Btr(v),  G\U'  =  r,  and  a(U,  T)  =  a(U',T). 

Proof.  Choose  r'  between  r  and  2 r  so  that  \n(Bri(v))\  <  A(3r.  Let  Uq  =  Bri(v),  and  G\Uq 
have  components  T^  . . . ,  Tp. 

Now  set  s  =  8a~1/2a~1/3A  and  ^  =  2s  +  2.  We  claim  that  all  but  one  of  the  components 
Tj  are  contained  in  B^r{v).  For  suppose  not;  then  for  i  ^  j  there  are  w  G  Tj  and  w'  G  Tj 
such  that  w  and  w'  are  each  at  distance  s  from  Uq,  Bs{w )  C  Tj,  and  Bs{w')  C  Tj.  Now 
consider  the  edge  cut  of  size  at  most  4/3 Ar  formed  by  5(Uo)]  one  of  the  two  spheres  Bs(w )  and 
Bs(w '),  say  the  latter,  is  contained  in  a  small  component  of  this  cut  in  G[B^r(v)}.  But  then 
the  semi-expansion  of  G[B^r(v)]  requires  that  4(3 Ar  >  ayJ\Bs(w')\,  which  is  a  contradiction 
since  by  Lemma  4.4  we  have  \Bs(w')\  >  as2. 

So  for  some  i,  only  T,;  is  not  contained  in  B^r{v).  Now  let  Tj, . . . ,  T'q  denote  the  compo¬ 
nents  of  G\U)  so  T,:  is  contained  in  one  of  these,  say  Tj,  and  T^, . . . ,  T^  are  all  contained  in 
B^r(v).  Thus  we  have 

U'  =  UU  |J  T'  C  B^r(v). 

j>  i 

In  particular,  U'  is  planar  since  f  r  <  A  log  n,  and  it  is  simple  since  G\U'  has  only  the 
component  Tj.  Thus  U'  satisfies  the  conditions  of  the  lemma.  ■ 

Finally,  we  show  a  general  property  of  planar  graphs  P[  with  small  face  size:  if  the 
distance  between  two  nodes  in  H  is  large,  than  the  value  of  any  edge  cut  which  contains 
both  in  the  same  segment  of  its  boundary  must  also  be  relatively  large. 
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Lemma  4.7  Let  H  be  a  planar  graph,  with  distinguished  faces  $1, . . . ,  $r  bowided  by  cycles 
Qi, ...  ,Qr  respectively.  Suppose  that  all  faces  other  than  $1, . . . ,  <3>r  are  bounded  by  at  most 
i  edges,  and  for  a  constant  d'  and  all  i  ^  j  we  have  d{Qi ,  Qj )  >  d' . 

Let  U  C  V  ( H )  and  v,  w  G  cr(U,  C )  for  some  component  C  of  G\U.  Then 

|5(i7)|  >  min (W1  d' ,  W1  d(v ,  w)) . 


Proof.  Let  S  =  cr(U,  C )  C  U.  In  the  graph  H[U\,  S  lies  on  a  single  facial  cycle  Q.  Traversing 
Q  in  a  clockwise  direction  starting  at  v,  we  encounter  faces  R\, . . . ,  Rp  whose  boundaries 
contain  vertices  both  of  U  and  of  H  \  U. 

Suppose  that  among  the  {/?,}  there  are  two  distinct  large  faces  &m  and  <Fm/.  Choose 
such  a  pair  for  which  Ra  =  <Fm,  Rb  =  and  Rc  {<hj}  for  a  <  c  <  b.  Let  P  denote  the 
corresponding  maximal  subpath  of  Q  whose  internal  vertices  are  incident  only  to  faces  Rc, 
for  a  <  c  <  b.  Then  among  every  l  consecutive  vertices  of  P,  there  must  be  one  incident 
to  an  edge  in  5(U);  since  \P\  >  d '  by  the  hypotheses  of  the  lemma,  this  implies  the  claimed 
bound. 

Otherwise,  there  is  a  single  large  face  <hm  among  the  {-R;};  note  that  may  appear 
several  times  on  the  traversal  of  Q.  Now  there  are  two  sub-paths  of  Q  from  v  to  w,  which 
we  denote  Po  and  P\.  Since  v,w  border  the  same  component  of  H  \  U,  the  face  <hm  does 
not  appear  in  a  traversal  of  one  of  P0  or  Pi  —  suppose  it  is  Pq.  So  as  above,  among  every 
t  consecutive  vertices  of  Po,  there  must  be  one  incident  to  an  edge  in  and  we  have 

\Po\>  d(v,  w).  m 

4.2  Related  Classes  of  Graphs 

In  this  section,  we  show  a  natural  construction  which  produces  uniformly  densely  embedded 
graphs;  it  is  an  extension  of  the  definition  of  geometrically  well-formed  graphs  in  our  earlier 
paper  [14],  The  material  in  this  section  is  independent  of  the  rest  of  the  paper. 

We  wish  to  define  a  notion  of  a  surface  being  locally  planar,  in  the  following  sense.  Let 
S  be  a  compact  orientablc  surface,  embedded  in  R3.  For  ifE,  let  B'd(x)  denote  the  set  of 
all  points  of  E  whose  distance  from  x  (as  measured  on  E)  is  at  most  d. 

Definition  4.8  A  set  X  C  E  is  (70,7i)-flat  for  some  positive  constants  7o,7i,  if  there  is  a 
E -disc  D  such  that 

(i)  X  C  D. 

(ii)  For  all  points  x  G  A"  and  s  >  0  such  that  B's{x )  C  X ,  the  surface  area  of  B's(x)  is  at 
least  70 s2  and  at  most  71  s2. 

(in)  For  all  E -disc  D'  such  that  D '  C  D,  if  D' ’ s  boundary  has  length  s,  then  the  surface 
area  of  D'  is  at  most  71  s2. 

We  say  that  E  is  (r,  70, 7i)-locally  flat  if  it  is  orientable,  and  for  all  x  G  E  the  set  B'r(x ) 

^  (70, 7i)-/^- 
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Of  course  all  these  properties  hold  if  E,  for  example,  is  the  unit  sphere  in  R3. 

Now  we  say  that  a  graph  is  locally  well-formed  if  it  is  drawn  on  a  locally  flat  surface,  and 
each  face  has  geometrically  about  the  same  (small)  size. 

Definition  4.9  A  graph  H  drawn  on  E  is  locally  well-formed  with  constant  parameters 
A,  £,  70, 71,  p0,  Pi  it  has  maximum  degree  A  and  there  is  an  r  >  0  such  that 

(i)  E  is  (r  log  n,  70, 71) -locally  flat, 

(ii)  The  maximum  number  of  edges  on  a  face  in  the  drawing  of  H  is  i,  and 
(in)  for  each  face  of  G  there  is  an  x  G  E  so  that  B'por(x)  C  C  B'pir(x). 

We  now  want  to  show  that  every  locally  well- formed  graph  is  uniformly  densely  embedded. 
To  show  this,  the  following  routine  lemma  is  useful:  in  Definition  4.1  it  is  enough  to  require 
semi-expansion  for  cuts  that  produce  only  two  components. 

Lemma  4.10  H  is  an  a-semi-expander  if  and  only  if  the  condition  of  Definition  f.l  holds 
for  all  sets  X  for  which  H[ X]  and  H  \  X  are  both  connected. 

Proof.  We  proceed  by  induction  on  the  number  of  connected  components  of  H[ X]  and  H\ X. 
Assume  H[ X]  is  not  connected,  and  let  Ti, . . .  ,TP  be  components.  Then  each  of  the  sets 
Ti, . . . ,  Tp  must  satisfy  Definition  4.1  by  the  induction  hypothesis.  From  this  we  get 

W*)l  =  E  wry  I  >  «E  v'ihi  >  “VTui  =  af\X\. 

i  i 

If  H  connected  but  H\ X  is  not,  and  each  connected  component  has  size  at  most  \\V (H) |, 
then  the  above  argument  applies  to  the  components  of  H\ X.  Otherwise,  apply  Definition  4.1 
to  the  cut  defined  by  the  single  large  component,  both  of  whose  sides  are  connected.  ■ 

Proposition  4.11  If  G  is  locally  well-formed  with  parameters  A,  £,  j0, 71,  p0,  pi,  then  there 
are  positive  constants  a  and  A  such  that  G  is  uniformly  densely  embedded  with  parameters 
a,  A,  A,  and  I. 

Proof.  Let  G  be  locally  well-formed  with  the  given  parameters.  Then  for  any  v  G  V, 
if  s  <  p^logn,  the  set  Bs(v)  is  contained  in  B'rlogn(v)  and  hence  in  a  E-disc.  Now  let 
X  C  Bs(v),  we  wish  to  show  that  it  satisfies  the  semi-expansion  requirement  in  G[Bs(v)]. 
By  Lemma  4.10,  we  may  assume  that  both  G[X\  and  G[Bs(v)  \  X]  are  connected.  Thus 
S(X)  lies  on  a  single  face  of  G[X].  Let  q  =  |5(A)|;  then  there  is  a  closed  curve  L  on  E  of 
length  at  most  p\rq  that  bounds  a  E-disc  containing  X.  Thus,  X  is  contained  in  a  disc  of 
area  at  most  71  p\r2qf.  But  each  face  in  G[X\  has  area  at  least  7 0p2r2,  so  X  has  at  most 
7iPi7o"Vo ' 2(12  feces,  and  hence  at  most  I  times  this  many  vertices.  ■ 

In  a  series  of  papers  proving,  among  other  things,  that  the  disjoint  paths  problem  for  a 
fixed  number  of  terminal  pairs  is  solvable  in  polynomial  time  [27],  Robertson  and  Seymour 
make  use  of  another  notion  of  “denseness”  of  surface  embeddings  —  namely  representativity. 
It  turns  out  that  our  definition  of  uniformly  densely  embedded  graphs  could  also  have  been 
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expressed  in  these  terms.  We  say  that  a  closed  curve  on  E  is  null-homotopic  if  it  is  homotopic 
to  a  point;  it  is  well-known  (see  e.g.  [25])  that  a  closed  curve  is  null-homotopic  if  and  only  if 
it  is  contained  in  a  E-disc.  Now  we  say  that  a  drawing  of  G  on  E  is  c-representative  [25,  26] 
if  any  non- null- homotopic  closed  curve  on  E  meets  the  drawing  of  G  at  least  c  times. 

In  this  terminology,  we  could  have  replaced  the  condition  that  each  G[Br{y)]  (r  <  A  log  n) 
be  contained  in  a  E-disc  by  the  condition  that  the  drawing  of  G  be  D  (log /^-representative. 
More  precisely, 

Proposition  4.12  If  G  satisfies  parts  (i)  and  (in)  of  Definition  f.2,  and  the  drawing  of 
G  is  (A  log  n) -representative,  then  there  is  a  constant  X'  such  that  G  is  uniformly  densely 
embedded  with  parameters  a,  A',  A,  and  I.  Conversely,  if  G  is  uniformly  densely  embedded 
with  parameters  a,  A,  A,  and  I,  then  there  is  a  constant  A'  such  that  the  drawing  of  G  is 
(A'  log  n) -representative. 

Proof.  The  converse  statement  is  easier.  If  G  is  uniformly  densely  embedded  with  parameters 
a,  A,  A,  and  I,  then  any  closed  curve  IZ  on  E  meeting  G  at  fewer  than  Alogn  vertices 
meets  it  only  at  vertices  contained  in  B\\ogn(v )  for  some  v  E  V.  Thus  IZ  is  contained  in  a 
E-disc  and  is  null-homotopic. 

Now  suppose  G  satisfies  parts  (i)  and  (iii)  of  Definition  4.2,  and  the  drawing  of  G  is 
(A  log /^-representative.  We  must  show  that  for  some  A',  every  G[B\ilognfv )]  is  drawn  in  a 
E-disc.  Choose  X'  <  (X  and  let  U  =  B\r\ogn(y)  for  some  v  E  V.  We  claim  that  every  simple 
cycle  of  G[U\  is  null-homotopic  in  E.  For  suppose  not,  and  choose  the  shortest  non-null- 
homotopic  cycle  Q  contained  in  G[U],  Say  for  simplicity  that  Q  contains  an  even  number  of 
vertices,  v0, . . . ,  v &, . . . ,  v2k  =  v0,  and  let  Q0  and  Qi  denote  the  two  sub-paths  of  Q  with  ends 
equal  to  vq  and  Vk-  Now  suppose  there  were  some  path  P  in  G[U\  with  ends  equal  to  vq  and 
Vk  of  length  less  than  k\  then  one  of  Qq  U  P  or  Q\  U  P  would  contain  a  non- null- homotopic 
simple  cycle  of  G[U\  of  length  less  than  2k,  contradicting  our  choice  of  Q.  Thus  k  <  2A'log  n, 
and  so  \Q\  <  4A'logn  <  Alogn.  Now  since  G  is  strongly  embedded,  there  is  a  simple  closed 
curve  IZ  on  E,  meeting  G  at  precisely  the  vertices  of  Q,  that  is  homotopic  to  Q  in  E;  but 
since  IZ  meets  G  fewer  than  Alogn  times,  it  is  null-homotopic  in  E.  This  contradicts  our 
assumption  that  Q  is  non-null-homotopic. 

Thus  G[U\  contains  only  null- homotopic  simple  cycles.  By  Theorems  (11.2)  and  (11.10) 
of  [25],  this  implies  that  G[U\  is  contained  in  a  E-disc.  ■ 


5  Densely  Embedded  Graphs:  The  Routing  Algorithms 

One  encounters  a  number  of  difficulties  in  extending  the  algorithms  of  Section  3  to  densely 
embedded  graphs  in  general.  Some  of  these  are  easily  taken  care  of —  for  example,  we  cannot 
define  “subsquares”  of  G  anymore;  but  we  can  use  balls  of  the  form  Br{y )  instead,  and  we 
have  seen  above  that  these  behave  in  much  the  same  way.  We  similarly  may  choose  a  maximal 
set  of  mutually  distant  balls  and  grow  enclosures  around  them.  The  major  problems  are  the 
following.  (1)  We  used  the  natural  crossbars  inside  a  mesh  for  routing;  do  these  enclosures 
have  similar  crossbars  inside  them?  (2)  Where  is  the  high-capacity  simulated  network  A/"? 
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To  build  crossbars  inside  the  enclosures  we  use  the  Okamura-Seymour  theorem  [20], 
analogously  to  a  construction  in  our  earlier  paper  [14].  To  define  the  high-capacity  simulated 
network  J\f  we  want  to  grow  the  enclosures  out  until  they  touch.  However,  at  this  point 
their  boundaries  might  not  be  “smooth”  enough  to  allow  us  to  build  crossbars  inside  them; 
additionally,  there  is  no  reason  for  enclosures  that  do  touch  to  have  fi(logn)  edges  in  their 
common  boundary. 

Nevertheless,  it  is  still  possible  to  build  a  simulated  network  A/",  as  follows.  We  grow 
enclosures  that  have  smooth  boundaries,  and  arc  large  enough  that  they  contain  large  cross¬ 
bars,  but  we  keep  them  mutually  distant  from  one  another.  Then  we  define  the  notion 
of  a  Voronoi  partition  of  G  to  allow  us  to  determine  which  clusters  are  “neighbors”;  we 
define  the  simulated  network  J\f  by  putting  fl(logn)  parallel  edges  between  neighboring  clus¬ 
ters  (whether  or  not  they  have  that  many  edges  in  the  common  boundary  of  their  Voronoi 
regions) . 

We  show  that  the  collection  of  these  parallel  edges  “represents”  the  graph  G  well  enough 
that  it  can  be  used  as  the  network  A/".  In  particular,  we  need  to  show  that  all  paths  accepted 
by  the  simulated  network  can  be  routed  in  G.  For  this  we  make  use  of  a  theorem  of  Schrijver 
[28];  we  show  that  there  exist  fl(logn)  paths  in  G  between  the  neighboring  enclosures,  such 
that  all  paths  between  all  pairs  are  mutually  disjoint. 

We  make  no  attempt  to  optimize  constants  here.  Set  Ao  =  A,  and  choose  positive 
constants  Ai,A2,  ...  so  that  AJ+i  <C  A  j  (the  exact  relationship  between  these  constants  is 
easy  to  determine  from  the  analysis  below).  A  connection  (sj,  A)  is  short  if  d(si,  U)  <  A2  log  n 
and  long  otherwise.  As  before,  we  handle  long  and  short  connections  separately;  for  now  we 
concentrate  on  pre-processing  the  graph  as  described  above  for  handling  long  connections. 

5.1  Building  the  Simulated  Network 

5.1.1  Clusters  and  Enclosures 

We  wish  to  choose  a  maximal  set  of  mutually  distant  vertices  around  which  to  grow  clusters. 
Let  Gr  denote  the  graph  obtained  from  G  by  joining  u  and  v  if  d(u,v )  <  r.  We  first  run 
Luby’s  randomized  maximal  independent  set  algorithm  [17]  in  CAslogn. 

Let  M  denote  the  resulting  Mis.  For  any  x  G  V,  some  vertex  within  A5logn  of  it  will 
enter  M  on  the  first  iteration  if  the  largest  number  chosen  in  -B2A3iogn(^)  is  chosen  by  a 
vertex  in  B\s\ogn(x).  This  happens  with  constant  probability,  by  Lemma  4.4.  Moreover,  if 
d(x,y )  >  A2  logn,  then  these  events  are  independent  for  x  and  y.  Thus, 

Lemma  5.1  Let  x,y  G  V  be  such  that  d(x,y )  >  A2  log  n.  Then  with  constant  probability 
there  are  u,v  G  M  such  that  d(x,u )  <  A5  log  n  and  d(y,v )  <  A5  log  n. 

Around  each  v  G  M  we  now  grow  a  cluster  of  radius  roughly  A5  log  n,  and  an  enclosure 
around  each  cluster,  with  “smooth”  boundaries.  We  need  the  following  facts.  Let  H  denote 
an  arbitrary  graph,  and  Q  a  simple  cycle  of  H .  For  u,  v  G  Q ,  let  g?q(m,  v)  denote  the  shortest- 
pat  h  distance  from  u  to  v  on  Q  —  that  is,  the  length  of  the  shorter  of  the  two  u-v  paths  on 

Q- 
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Figure  2:  Building  the  simulated  network 

Definition  5.2  We  say  that  Q  is  e-smooth  if  for  all  u,v  E  Q  we  have  £cIq(u,v )  <  d(u,v). 

Definition  5.3  If  U  and  W  are  two  subsets  ofV(H),  we  say  that  U  is  e' -close  to  W  if  for 
each  ueU  there  is  a  w  E  W  such  that  d(u,w )  <  e'|bF|. 

The  following  fact  is  quite  similar  to,  but  more  general  than,  Theorem  4.4  of  onr  earlier 
paper  [14];  the  proof  is  very  similar  as  well.  In  effect,  it  says  that  given  a  cycle  Q  in  a 
planar  graph  H  that  encloses  (in  the  sense  of  homotopy)  the  “hole”  formed  by  some  internal 
face,  then  for  a  small  e  >  0  we  can  find  a  cycle  Q'  of  no  greater  length  that  is  e-close  to  Q, 
n(e)-smooth,  and  also  encloses  this  hole.  See  Figure  3. 

We  will  use  this  theorem  to  smooth  out  the  boundaries  of  the  clusters  and  the  enclosures 
around  them. 

Theorem  5.4  For  each  £  >  0  the  following  holds.  Let  Si  be  a  compact  surface  (possibly 
with  boundary),  H  a  graph  embedded  on  S1?  and  Q  a  simple  cycle  of  H  that  is  non-null- 
homotopic  on  Si.  Then  in  polynomial  time  one  can  find  an  -^—-smooth  simple  cycle  Q'  such 
that 
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Figure  3:  Smoothing  a  cycle 


(i)  \Q'\  <  \Q\, 

(ii)  Q'  is  E-close  to  Q,  and 

(in)  Q'  is  also  non-null-homotopic  on  E^ 

Proof.  For  u,v  G  Q,  let  [u,v\q  denote  the  shorter  of  the  two  u-v  paths  contained  in  Q  (ties 
broken  arbitrarily),  and  let  e  =  If  Q  is  not  e-smooth,  then  there  are  u,v  G  Q  such  that 

edQ(u,v )  >  d(u,v).  (3) 

Moreover,  we  can  efficiently  ford  such  a  u  and  v  so  that  there  is  a  shortest  u-v  path  Puv  in 
PI  that  is  vertex-disjoint  from  Q  (for  example,  the  pair  u,v  satisfying  (3)  for  which  \PUV\  is 
minimum) . 

Now  one  of  the  two  simple  cycles  [u,  v\qUPuv  and  (Q\  [u,  v]q)UPuv  is  not  null-homotopic 
on  Ei;  and  each  is  shorter  than  Q.  We  thus  update  Q,  replacing  it  with  the  cycle  from  among 
these  two  that  is  not  null-homotopic. 

We  now  iterate  this  process  of  “slicing  off”  parts  of  Q  using  short  paths  through  H .  Since 
the  length  of  the  cycle  decreases  with  each  iteration,  this  process  must  terminate  in  a  cycle 
Q'  that  is  e-smooth.  Moreover,  each  iteration  maintains  the  invariant  that  the  current  cycle 
is  non- null-homotopic  on  Ei.  Thus,  we  only  have  to  verify  that  the  final  cycle  is  e-close  to 

Q- 
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This  is  clearly  true  after  the  first  iteration:  since  \PUV\  <  £cIq(u,v )  <  e\Q\,  every  vertex 
on  the  updated  cycle  can  reach  a  vertex  of  Q  by  a  path  of  length  at  most  e\Q\.  Now,  let  Qi 
denote  the  cycle  obtained  after  i  iterations  of  slicing  off.  As  long  as  portions  of  Q  remain 
on  Qi,  we  say  that  we  are  in  the  “first  phase”;  other  phases  will  be  defined  below.  In  the 
first  phase,  Qi  consists  of  alternating  intervals  Qa,  Pr\ ,  Qi 2, . . . ,  Qir,  Pir,  where  QtJ  C  Q,  and 
the  interval  Q'V]  of  Q  lying  between  QVJ  and  Qi,j+i  has  been  sliced  off  by  P^.  We  can  show 

by  induction  on  the  number  of  iterations  that  PtJ  <  e  Q[j  —  as  was  true  after  the  first 
iteration. 

This  is  done  by  the  following  case  analysis.  In  the  (i  +  l)st  iteration,  we  hnd  a  new  path; 
there  are  three  cases  to  consider. 

1.  One  end  of  P  lies  on  Qij  and  the  other  on  Qa, ,  where  possibly  j  =  k.  Then  the 
property  clearly  continues  to  hold,  since  |P|  is  at  most  £  times  the  number  of  current 
cycle  vertices  cut  off,  which  is  in  turn  at  most  the  number  of  original  vertices  of  Q 
between  the  endpoints  of  P. 

2.  One  end  of  P  lies  on  Pt3  and  the  other  on  Qik  (so  Pl}  is  lengthened).  Suppose  that  the 
amount  of  original  cycle  cut  off  in  addition  to  Qt-  is  equal  to  x,  and  the  amount  of  Pi3 


that  is  cut  off  by  P  is  y.  Then  if  Pl+l  j 

denotes  PV]  after  this  iteration,  we  have 

m 

< 

E  Qij 

\p\ 

< 

£(x  +  y) 

\Pi+l,j\ 

= 

\Pij  \  +  \p  \  —  y 

<  e(| Q,ij\+x  +  y)-y 

<  £(\Q'ij\+x) 

3.  One  end  of  P  lies  on  PtJ  and  the  other  lies  on  Pik  (so  P  glues  some  of  the  new  paths 
together).  There  are  two  subcases. 

(i)  j  =  k.  Then  Pt]  goes  down  while  \QZ3  \  is  not  affected,  so  the  property  still  holds. 

(ii)  j  7^  k.  Again  suppose  that  the  amount  of  original  cycle  cut  off  in  addition  to  Qh 
and  Q'ik  is  equal  to  x,  the  amount  of  PtJ  cut  off  by  P  is  y,  the  amount  of  Pik  cut 
off  by  P  is  z,  and  the  new  interval  is  denoted  Pi+ij-  Then 


\Pn\ 

< 

eQij 

P^ 

< 

z\Q'ik\ 

\p\ 

< 

E(x  +  y  +  z) 

Pi+l,j\ 

= 

|  Pij  |  +  \P\  +  \Pik\  ~  V  —  Z 

< 

z{Q'ij  +x  +  y  +  z+  \Q'ik\)  —  y  —  z 

< 

£{Q'ij  +x  +  \Q'ik\ ) 
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If  the  iterations  come  to  an  end  before  the  end  of  the  first  phase,  then  indeed  Q'  is  e-close 
to  Q  —  any  vertex  on  P \j  can  reach  Q  by  a  path  of  length  at  most  e  Q't]  <  e\Q\.  Otherwise, 
consider  the  iteration  in  which  the  first  phase  comes  to  an  end.  By  analogous  arguments, 
we  obtain  a  cycle  Q 1  such  that  \Ql\  <  s\Q\  and  every  vertex  on  Q 1  can  reach  Q  by  a  path 
of  length  at  most  s\Q\. 

Each  phase  now  proceeds  exactly  like  the  previous  one,  except  that  it  begins  with  a  cycle 
whose  length  has  been  reduced  by  at  least  a  factor  of  e.  Thus  when  the  process  terminates, 
all  vertices  on  Q'  will  be  able  to  reach  Q  by  a  path  of  length  at  most 

OO 

1=1 

Thus  Q'  is  e-close  to  Q.  ■ 


We  use  the  following  procedure  to  build  the  clusters  and  the  enclosures  around  each  node 
v  in  M.  Let  Kv  =  B\5iogn(v). 

(i)  Choose  a  radius  r  between  2A3  log  n  and  3A3  log n  so  that  |7r(Sr(u))|  <  9/3A.5  logn.  Set 
Cv  =  Br{y). 

(ii)  Now  extend  Cv  to  a  simple  set  as  in  Lemma  4.6;  since  A3  >  2cA5,  no  Cv  is  grown 
enough  that  it  overlaps  any  other  by  this  process. 

(iii)  We  now  apply  the  ^-smoothing  algorithm  of  Theorem  5.4  to  the  facial  cycle  Qv  of 
G[CV\  containing  tt(Cv).  Here  Hv  =  G[B\3  iogn(u)]  \  K°  plays  the  role  of  H,  and  the  cylinder 
formed  by  removing  the  portions  of  E  on  which  G[K°\  and  G  \  B\3\ogn(v)  are  drawn  plays 
the  role  of  E3.  Now  for  a  constant  £  the  resulting  cycle  Q'v  is  e-smooth  in  this  subgraph  Hv 


of  G,  and  it  is  also  ^ly-close  to  Qv.  If  we  choose  e  <  18^+1 ,  then  since  Qv  is  initially  at  least 
7^j\Qv\  away  from  Kv ,  we  know  that  any  path  with  both  ends  on  Q'v  that  passes  through  Kv 


must  have  length  at  least  ^\QV\  >  ^g\Q'v\-  Thus,  there  are  no  “short  cuts”  between  vertices 
of  Q'v  that  make  use  of  Kv ;  hence  Qv  is  in  fact  e-smooth  in  G. 


The  smooth  cycle  Q'v  encloses  a  set  S  of  vertices  containing  Kv.  Update  Cv  to  be  this 


set  S. 


We  now  grow  an  enclosure  Dv  D  Cv  by  the  same  three-step  process,  except  that  we  now 
use  the  constant  A4  in  place  of  A5,  and  the  set  C°  in  place  of  K°  Thus,  we  have  clusters 
of  radius  ~  Aslogn,  enclosures  of  radius  ~  A4log  n,  and  they  are  separated  by  a  mutual 
distance  of  rs  A3  log  n. 

Following  the  algorithm  of  Section  3,  we  now  must  build  crossbar  structures  in  the  en¬ 
closures  to  replace  the  natural  crossbars  of  the  mesh.  We  build  crossbars  using  an  extension 
of  the  OkamuraUSeymour  theorem  [20]  due  to  Frank  [9],  along  the  same  lines  as  was  done 
in  [14],  To  be  precise, 


Definition  5.5  If  X  C  V ,  we  say  a  crossbar  anchored  in  X  is  a  set  of  edge-disjoint  paths, 
each  with  at  least  one  end  in  X,  such  that  every  pair  of  paths  meets  in  at  least  one  vertex. 


Let  av  =  n(Cv),  and  recall  that  by  Q'v  we  mean  the  facial  cycle  of  G[CV\  containing  av . 
Analogously,  let  rv  =  7 r(Dv),  and  Lpv  the  facial  cycle  of  G[DV\  containing  tv .  We  wish  to 
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build  a  crossbar  in  G[DV  \  C°],  anchored  in  av  U  rv,  of  size  at  least  a  constant  fraction  of 
(7V  U  tv | .  For  a  large  enough  constant  k  depending  on  e,  we  choose  a  set  a'v  ot\av\/ k  vertices 
on  <jv  spaced  about  k  apart,  and  a  set  t'v  of  \tv\/k  on  r„  spaced  about  k  apart. 

Lemma  5.6  There  is  a  crossbar  anchored  in  a'v  U  t'v,  such  that  each  vertex  of  a'v  U  r'v  is  the 
endpoint  of  a  distinct  path  of  the  crossbar. 

Proof.  Consider  the  planar  graph  G[DV  \  C°].  This  graph  has  only  two  large  faces  —  the 
outer  face  bounded  by  (pv,  and  the  inner  face  left  by  the  deletion  of  C°.  We  find  a  shortest 
path  P*  in  the  planar  dual  graph  whose  endpoints  are  equal  to  these  two  large  faces,  and 
we  delete  the  edges  used  by  this  path.  Denote  the  resulting  graph  by  Gv ;  note  that  it  has 
only  a  single  large  face,  bounded  by  a  cycle  ip'v  which  contains  a'v  U  r'v. 

We  claim  that  the  cycle  ip'v  is  e'-smooth  in  the  graph  Gv ,  where 

s'  =  min(^£,  (27 

To  see  this,  suppose  that  P  is  a  path  in  Gv  with  endpoints  u,w  G  ip'v.  If  u  and  w  both  belong 
to  Q'v ,  or  both  belong  to  tpv,  then  this  follows  from  the  e-smoothness  of  these  two  cycles  in 
G[DV  \  C°\.  If  one  belongs  to  each,  then  |P|  >  (A4  —  4A5)  logn,  while  \<p'v\  <  9/3£(2A4  +  A5), 
and  again  the  bound  follows.  Finally,  if  both  u  and  w  lie  on  the  short  dual  path  P*, 
then  in  fact  \P\  >  dtp'v(u,w),  while  if  exactly  one  lies  on  P*  then  it  is  easily  verified  that 
1^1  >  | ed^faw). 

Write  X  —  a'v  U  r'.  Let  us  assume  for  simplicity  that  |A"|  is  odd.  Now  for  u  G  X,  define 
u+  to  be  the  vertex  on  ip'v  that  is  | k,  steps  clockwise  from  u,  and  write  X+  =  {u+  :  u  G  X}. 
Now  \X  U  A"+|  is  even,  so  we  can  pair  each  «GlU  A"+  with  its  unique  “antipodal”  point 
u  e  XU  X+  under  the  clockwise  ordering  of  ip'v.  Note  that  vertices  in  X  are  paired  with 
vertices  in  A"+,  and  vice  versa.  We  now  define  a  disjoint  paths  problem  in  Gv ,  with  the  set 
of  terminal  pairs  %  equal  to  {(u,u)  :  u  G  X}.  Note  that  all  terminals  are  at  least  from 
one  another  on  the  cycle  <p'v. 

We  now  want  to  show  that  Tv  is  realizable  in  Gv.  First,  say  that  a  cut  is  non-trivial  if 
it  separates  at  least  one  pair  of  terminals.  We  are  dealing  with  a  disjoint  paths  problem 
in  a  planar  graph  with  all  terminals  on  the  outer  face.  Moreover,  every  node  of  Gv  not  on 
the  outer  face  has  even  degree.  In  such  a  case,  Frank’s  extension  of  the  Okamura-Seymour 
theorem  [9,  20]  says  that  following  strict  cut  condition  is  sufficient  for  the  realizability  of  Tv 
-  every  non-trivial  cut  has  more  capacity  than  the  number  of  terminal  pairs  it  separates. 

It  is  enough  to  consider  non-trivial  cuts  of  the  form  S(U)  with  both  GV[U]  and  GV\U 
connected.  For  such  a  set  U,  there  must  be  two  vertices  v,w  G  7r (U)  such  that  v,w  G  ip'v. 
Suppose  that  the  distance  from  v  to  w  in  Gv  is  d;  then  by  Lemma  4.7,  we  have  |5(f7)|  > 

Since  the  facial  cycle  ipv  is  A-smooth,  the  number  of  terminal  pairs  disconnected  by  S(U) 
is  at  most  2 e'^d/n.  Thus  taking  k  >  2G'~l  ensures  that  the  strict  cut  condition  will  be 
satisfied. 

Finally,  observe  that  the  edge-disjoint  paths  in  a  realization  of  %,  provide  the  crossbar 
required  by  the  lemma,  since  each  pair  of  paths  must  meet  at  some  vertex  of  Gv .  m 
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In  the  crossbar  just  constructed,  let  Y™  denote  the  path  with  an  endpoint  equal  to  u, 
where  u  G  a'v  U  t'. 

5.1.2  The  simulated  network  M 

We  now  construct  the  simulated  network ;  the  nodes  of  this  network  are  the  clusters,  which  we 
represent  by  the  vertices  in  M .  We  define  a  neighbor  relation  on  the  clusters  using  the  notion 
of  a  Voronoi  partition ;  two  clusters  will  be  joined  by  an  edge  in  J\f  if  they  are  neighbors  in 
this  sense. 

Let  H  be  a  graph  and  S  C  V(H).  We  fix  a  lexicographic  ordering  ■<  on  the  elements  of 
S.  For  s  G  S.  let 

Us  —  {v  G  G  :  Vs'  G  S  :  d(v,  s )  <  d(v,  s')  and  Vs'  V  s  :  d(v,  s )  <  d(v,  s')}. 

That  is,  Us  is  the  set  of  vertices  that  are  at  least  as  close  to  s  as  to  any  other  s',  with  ties 
broken  based  on  V 

Definition  5.7  The  Voronoi  partition  V(H,  S )  of  H  with  respect  to  S  is  the  partition  {Us  : 
seS}. 

The  following  fact  is  immediate. 

Lemma  5.8  For  each  s  G  S ,  H[US]  is  connected. 

Proof.  Suppose  v  G  U.s\  we  claim  that  any  shortest  s-v  path  P  is  contained  in  Us.  For 
suppose  not,  and  let  v'  G  P  be  the  closest  vertex  to  s  that  lies  in  Us'  for  some  s'  7^  s. 
Then  d(s',v )  <  d(s,v),  and  in  fact  d(s',v )  <  d(s,v )  if  s  V  s'.  It  follows  that  v  G  Us ?,  a 
contradiction.  ■ 

We  can  now  build  a  graph  A f(H,  S )  on  the  vertices  in  S,  joining  two  if  their  Voronoi  cells 
share  an  edge. 

Definition  5.9  The  neighborhood  graph  of  S  in  H,  denoted  A f(H,S),  is  the  graph  with 
vertex  set  S,  and  an  edge  (s,  s')  iff  there  is  an  edge  of  H  with  endpoints  in  Us  and  Usi . 

The  simulated  graph  we  use  will  be  the  neighborhood  graph  A f(G,  M )  with  every  edge 
given  capacity  ~  A6  log  n.  Let  V  and  J\f  denote  V(G,M)  and  A f(G,M)  respectively,  and 
A/"(7)  the  graph  A/"  in  which  each  edge  is  given  capacity  7. 

The  following  two  facts  about  A/"  are  easy  to  establish.  First,  by  the  maximality  of  M, 
we  have 

Lemma  5.10  For  all  v  G  M,  Uv  C  B\3iogn(v). 

Proof.  Suppose  u  G  Uv  but  d(v,u )  >  A3  log n.  Then  d(v',u )  >  A3  log  n  for  all  v'  G  M;  this 
contradicts  the  fact  that  M  is  a  maximal  independent  set  in  (7Aslogn.  ■ 

This,  along  with  Lemma  4.4,  implies 
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Lemma  5.11  The  degree  of  a  vertex  in  J\f  is  at  most  A'  <  16a  1(3. 

Proof.  Let  U  denote  the  neighbors  of  v  in  A/",  including  v  itself.  Then  by  Lemma  5.10, 

U  uu  P-  -^2A3log  n(v)i 

u&U 


and  hence 

^\Uu\<  4/7A2  log2  n. 

uGU 

But  around  each  u  G  U  there  is  a  disjoint  ball  of  radius  |A3logn,  which  contains  at  least 
dA|log2  n/4  vertices.  Thus  \U\  <  16d_1/3.  ■ 

5.1.3  Inter-Cluster  Paths 

The  goal  of  this  part  is,  for  a  constant  A6,  to  construct  A6  log  n  disjoint  paths  between  each 
pair  of  enclosures  DV,DW  where  (v,w)  is  an  edge  in  A f.  This  will  allow  us  to  convert  a 
routing  in  the  simulated  network  log  n)  into  actual  disjoint  paths  in  G.  Recall  that  the 
outer  facial  cycle  of  G[DV\  is  denoted  ipv,  and  it  contains  a  set  r'v  of  vertices  evenly  spaced 
at  distance  k. 

Theorem  5.12  There  exist  vertex- disjoint  paths  in  G,  each  with  ends  in  sets  r'  mid  oth¬ 
erwise  disjoint  from  all  Dv,  such  that  for  (v,w)  G  A f,  there  are  at  least  A&  log  n  such  paths 
with  one  end  in  r'  and  the  other  in  t'w. 

Proof.  The  proof  is  based  on  the  following  theorem  of  Schrijver  [28].  Let  E3  be  a  surface 
(possibly  with  boundary),  H  a  graph  embedded  on  E1}  and  { A 1  :  i  —  1, . . . ,  k}  a  set  of 
disjoint  curves  on  Ei  each  of  which  is  either  closed  or  anchored  at  vertices  of  PI  on  the 
boundary  of  Ei.  The  problem  is  to  find  vertex-disjoint  paths  and  cycles  {P1}  in  G  so  that 
Pl  is  homotopic  to  A 1  for  each  i. 

Call  a  collection  of  curves  7 Z  on  Ex  essential  if  it  consists  either  of  a  single  closed  curve 
that  is  not  null-homotopic,  or  it  is  a  finite  union  of  curves  each  with  endpoints  on  the 
boundary  of  Ei.  Schrijver  [28]  proves  that  such  vertex-disjoint  paths  and  cycles  exist  if  for 
each  essential  collection  of  curves  7 Z,  there  are  curves  {B1},  where  Bl  is  homotopic  to  A ®, 
such  that  1Z  intersects  the  drawing  of  H  more  than  it  intersects  the  curves  { B *}.  (The  main 
result  of  [28]  is  in  fact  a  necessary  and  sufficient  condition;  this  weaker  statement  suffices  for 
our  purposes.  Additionally,  [28]  is  stated  for  the  special  case  of  surfaces  without  boundary, 
but  the  extension  we  use  here  follows  immediately  from  [28].) 

Say  that  a  curve  is  G'-normal  if  it  meets  the  drawing  of  G  only  at  vertices,  and  define  its 
G'-length  to  be  the  number  of  times  it  meets  the  drawing.  For  each  v,  v'  that  are  neighbors 
in  TV",  we  draw  a  G'-normal  arc  Avv>  on  E  with  endpoints  v  and  v' .  We  can  ensure  that  all 
these  arcs  are  disjoint,  since  each  Uv  is  connected,  and  for  each  (v,v')  G  E(Af),  there  is  at 
least  one  edge  of  G  with  endpoints  in  Uv  and  Uvt.  Choose  a  small  constant  A6  <  \r'v |/A'  (say, 
less  than  |a2dA'~2A_2A3  A2;  the  reason  for  this  will  become  clear  below).  Now  suppose  we 
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Figure  4:  Pulling  out  a  crossing 

have  in  fact  log  n  copies  of  each  Avv',  all  running  “parallel”  to  one  another.  By  pushing 
them  apart  appropriately,  we  can  assume  that  each  arc  runs  through  a  different  vertex  in  t'v 
and  t',.  Let  Auut  denote  the  G'-normal  arc  that  runs  through  u  G  r'v  and  v!  G  t'v,. 

Define  G'  to  be  the  graph  obtained  by  deleting  D°  for  each  v  E  M.  We  now  cut  E 
along  the  facial  cycles  (pv  to  obtain  S',  a  surface  with  boundary.  Note  that  G'  is  properly 
embedded  on  S'.  The  theorem  is  now  a  consequence  of  the  following  claim.  ■ 


Claim  5.13  There  exist  vertex- disjoint  paths  Puu ,  in  G'  such  that  Puui  is  homotopic  to  Auu>. 


Proof.  If  curves  TZ  and  7 Z1  are  homotopic,  we  write  TZ  ~  TZ' .  We  extend  this  notation  to 
finite  collections  of  curves  TZ ,  TZ  in  the  obvious  way.  Following  notation  of  Schrijver  [28], 
the  number  of  crossings  of  TZq  and  TZi  is  denoted  cy(JZo,TZi),  and  we  define 

miner (7Z0,  TZ\ )  =  mm.{cr{JZ!0,TZ!1)  :  IZq  ~  TZ'0,TZi  ~  TZ'i). 

By  the  theorem  of  Schrijver  [28]  given  above,  the  desired  paths  exist  if  for  each  essential 
collection  of  curves  TZ  on  S',  one  has 

ct(JZ,G')>  ^  miner (TZ,Auur).  (4) 

(' U,U‘ ') 


Note  that  in  verifying  this  inequality,  we  may  assume  TZ  is  G'-normal,  and  that  TZ  has  no 
self-crossings. 
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For  a  collection  of  curves  TZ,  define  its  index  to  be 

cr(7 Z,G')~  Y  miner (TZ,AUU>). 

(u,u') 


So  it  is  enough  to  consider  collections  of  curves  TZ  whose  index  is  minimum  in  their  homotopy 
class,  and  to  show  that  such  TZ  in  fact  have  positive  index. 

Set  r  =  !A~1ad1/2A4  log  n.  We  claim  that  no  G- normal  closed  curve  of  (S'- length  less  than 
2 r  can  enclose  a  set  of  the  form  Dv.  For  if  it  did,  then  Dv  could  be  disconnected  by  an  edge 
cut  of  size  less  than  acA^2\\  logn,  which  is  not  possible  since  \DV\  >  dA2  log2  n.  From  this 
it  follows  that  any  G'-normal  closed  curve  of  G'-length  less  than  2 r  must  be  null-homotopic 
in  S'. 

Note  that  we  can  view  the  expression  cr(7 Z,G')  —  J2(u,u')  cr(7\L  Auu/)  as  a  sum  over  the 
finitely  many  arc-components  of  TZ;  if  the  value  of  this  expression  is  not  positive,  we  show  how 
modify  the  curves  Auu/  so  that  it  increases.  We  do  this  by  considering  each  arc-component 
of  TZ  in  turn.  Let  TZ  denote  a  single  arc-component  of  TZ;  we  consider  two  cases,  based  on 
the  G'-length  of  TZ. 

Case  1.  ct(1Z,G')  <  r.  Then  TZ  must  have  both  endpoints  on  the  same  facial  cycle 
(it  is  too  short  to  touch  two  such  cycles,  and  if  it  were  a  closed  curve  it  would  have  to  be 
null-homotopic,  by  the  above  argument.)  But  then  it  is  easy  to  produce  arcs  {A'uu}  for  which 
cv(TZ,  G')  >  J2(u,u')  cr(7 Z,  A'uu,)  since  <pv  is  e-smooth. 

Case  2.  cr(JZ,  G")  >  r.  Again,  we  just  have  to  exhibit  arcs  A'uu,  ~  Auui  lying  on  S'  so 
that 

a(R,G')>  Y.  cr(K,4g,),  (5) 

(u,u') 

without  increasing  the  number  of  crossings  of  these  arcs  with  the  other  components  of  TZ.  If 
the  set  {Auui}  satisfies  (5),  we  are  done;  otherwise,  we  show  how  to  modify  this  set  of  arcs 
so  that  it  does.  See  Figure  4. 

If  the  set  {Auui}  does  not  satisfy  Inequality  (5),  then  there  is  some  interval  TZ'  of  TZ  of 
G'-length  r  for  which  (5)  is  violated.  Let  us  consider  such  an  TZ' . 

Observe  that  each  arc  Auut  has  G'-length  at  most  2A3  log  n,  and  hence  at  most  A'2A6  log  n 
of  these  arcs  can  meet  TZ' ,  since  at  most  A'2  pairs  of  clusters  have  at  least  one  end  close 
enough  to  TZ' .  Now  suppose  the  total  number  of  crossings  of  these  arcs  with  TZ'  exceeds 

(2A:i  log  n) ( A'W  log  re)  <  ct(r,  g) 
r 

Then  some  arc  Auu/  meets  TZ'  more  than  2 A3  log  n/r  times,  and  hence  the  interval  of  Auui 
between  some  pair  of  consecutive  crossings  with  TZ'  has  G'-length  less  than  r. 

Suppose  that  this  pair  of  consecutive  crossings  occurs  at  vertices  w  and  w' .  Let  TZ"  denote 
the  G'-normal  curve  formed  from  this  interval  of  Auui  and  the  portion  of  TZ'  between  w  and 
w'.  TZ"  has  G'-length  less  than  2 r,  and  so  it  must  be  null- homotopic  by  the  argument  given 
above. 
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Now,  since  1Z  has  minimum  index  over  all  curves  in  its  homotopy  class,  the  portion  of 
Auui  between  w  and  w'  meets  G'  at  least  as  many  times  as  the  portion  of  TZ'  between  w  and 
w' .  We  can  therefore  modify  Auui  so  that  it  runs  along  TZ'  for  this  interval.  This  does  not 
increase  the  G'-lcngth  of  Auui ;  and  it  decreases  the  number  of  crossings  of  TZ  —  as  well  as 
TZ  (since  it  has  no  self-crossings)  —  with  Auu/. 

Thus  this  process  terminates;  when  it  does,  we  have  a  set  of  arcs  {A'uu,}  for  which 
Inequality  5  holds.  ■ 

Let  us  denote  one  such  path  with  ends  u  and  v!  by  Zuut.  Moreover,  we  have  u  G  r'v 
and  v!  G  t'v ,,  and  they  are  ends  of  paths  Yff  and  Yf)  respectively.  Denote  by  Zuu>  the 
concatenation  of  the  three  paths  Yff,  Zuu>,  and  Yft  . 

5.2  The  On-Line  Algorithm 

With  the  simulated  network  A/"  in  place,  the  routing  algorithms  themselves  are  essentially 
the  same  as  those  for  the  mesh;  here  we  just  describe  what  must  be  modified. 

First  of  all,  the  analogue  of  Lemma  3.10  is  the  following. 

Lemma  5.14  Let  r  <  Ailogn,  U  C  Br(v)  for  some  v  G  V,  and  T  a  set  of  terminal  pairs 
in  U.  Then  the  maximum  size  of  a  subset  of  T  that  is  realizable  in  B8^r(v)  is  within  a 
constant  factor  of  the  maximum  size  of  a  subset  of  T  that  is  realizable  in  G. 

Proof.  First  choose  a  radius  r'  between  2 r  and  3r  for  which  \n(Br(v)\  <  9 f3r' .  Then 
construct  a  simple  set  extension  of  Bri(v )  as  in  Lemma  4.6;  and  e-smooth  its  outer  facial 
cycle  to  obtain  a  set  U'  D  U  contained  in  B^r{v).  Let  U"  C  B8^2r  denote  a  simple  set 
extension  of  B8^r(v) ,  as  in  Lemma  4.6.  For  a  constant  we  can  pick  a  set  S  of  vertices  on 
the  outer  facial  cycle  of  U'  spaced  k'  apart,  and  use  Frank’s  theorem  [9]  as  in  Lemma  5.6 
to  construct  a  set  of  edge-disjoint  paths  connecting  “antipodal”  pairs  in  S,  such  that  all 
paths  stay  within  U"  \  U' .  Note  that  we  must  take  care  to  ensure  that  the  parity  condition 
is  met,  since  the  outer  facial  cycle  of  U"  can  contain  odd-degree  vertices.  To  do  this  we 
remove  sub-paths  of  this  cycle  between  consecutive  pairs  of  the  (necessarily  even  number  of) 
odd-degree  vertices;  U"  is  large  enough  that  the  strict  cut  condition  will  remain  satisfied. 

Consider  the  set  of  paths  in  a  realization  of  a  maximum-size  subset  of  T'  in  G.  Of  the 
paths  that  meet  we  can  select  a  constant  fraction  of  paths  such  that  the  the  pairs  of 

first  and  last  intersections  with  n (U')  can  be  connected  along  the  outer  facial  cycle  of  U’  to 
different  vertices  in  S.  We  can  then  use  the  crossbar  of  the  previous  paragraph  to  connect  all 
of  these  pairs  together;  the  resulting  paths  are  within  a  constant  fraction  of  the  maximum 
number  achievable  in  G,  and  they  do  not  leave  U".  ■ 

The  algorithm  for  short  connections  is  now  as  follows.  We  run  a  randomized  version  of 
Luby’s  algorithm,  this  time  in  GAllogn.  Let  M'  denote  the  resulting  Mis.  With  constant 
probability,  both  ends  of  a  short  connection  are  within  ^2 Ailogn  of  the  same  v  G  M',  as 
in  Lemmas  3.11  and  5.1.  We  now  let  Uv  denote  BXl  iogn/i6$2('i0  and  only  route  connections 
both  of  whose  ends  lie  in  the  same  Uv.  To  route  such  connections,  we  run  the  algorithm  of 
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Proposition  3.9  in  each  BiXilogn(v){  by  Lemma  5.14,  this  is  within  0(log n)  of  optimal  in 
each  Uv. 

For  long  connections,  we  run  the  same  simulation  as  in  the  case  of  the  mesh,  this  time 
in  the  graph  A/^Aelogn).  The  analogue  of  Lemma  3.6  holds  exactly  as  before,  as  does 
Lemma  3.7  —  making  use  of  the  paths  Zuu/,  since  all  such  paths  incident  to  the  same 
enclosure  must  cross.  Finally,  the  analogue  of  Lemma  3.8  is  the  following. 

Lemma  5.15  There  is  a  constant  7  such  that  for  any  realizable  subset  T  of  T,  can 

be  routed  in  7  log  n) . 

Proof.  Set  7  =  9/3(A/3A5  +  A3).  For  each  path  P  in  the  optimal  routing,  construct  the 
following  path  for  ( u^Vi )  in  J\f  —  when  P  crosses  from  Uw  into  IAW> ,  add  an  edge  from  w 
to  w' .  Now  consider  how  many  paths  in  our  constructed  routing  use  the  edge  ( w,w ').  Each 
such  corresponds  to  a  path  in  the  original  routing  that  used  an  edge  in  5{UW)UW').  We  can’t 
bound  the  size  of  this  set  directly,  since  it  could  be  quite  “meandering.”  But  consider  the 
following  argument.  5{UW)UW')  C  -B2A3iogn(r)  for  some  vertex  x  G  G;  thus  there  is  some  r 
between  2A3logn  and  3Aslogn  so  that  |7r(i?r(a;))|  <  9/9A3logn.  So  at  most  this  many  paths 
with  both  ends  more  than  r  away  from  x  can  use  edge  in  5{UW)UW').  But  closer  than  this, 
there  are  at  most  A'3  clusters,  each  of  which  is  the  origin  of  at  most  Qf3\^\ogn  paths.  ■ 

Thus  the  optimum  in  G  is  bounded  by  the  fractional  optimum  in  A/"(7logn),  which  is  at 
most  a  constant  factor  more  than  the  fractional  optimum  in  Af(\a  logn),  which  is  at  most  an 
O(logn)  factor  more  than  the  number  of  pairs  routed  by  the  on-line  algorithm  in  G.  Thus 
we  have 

Theorem  5.16  The  on-line  algorithm  is  O (logn) -competitive  in  any  uniformly  densely  em¬ 
bedded  Eulerian  graph  G. 

5.3  The  Off-Line  Algorithm 

The  off-line  algorithm  too  is  essentially  the  same,  now  working  with  the  larger  simulated 
network  logn).  The  only  change  required  is  in  the  proof  of  Lemma  3.15.  Here,  we  are 

no  longer  able  to  talk  about  “overfull  rectangles” ;  however,  we  can  define  a  round  cut  to  be 
a  set  of  the  form  Br(w)  fl  Cv.  Since  a  given  u  G  Cv  is  only  contained  in  0(r2)  round  cuts  of 
radius  r,  the  following  lemma  establishes  that  round  cuts  can  be  used  instead  of  rectangles 
in  the  analogue  of  Lemma  3.15  for  densely  embedded  graphs. 

Lemma  5.17  Let  G'v  denote  G[CV]  with  an  additional  vertex  zv  joined  by  an  edge  to  each 
vertex  in  a'v.  Then  there  is  a  constant  fi  such  that  for  every  U  C  G'v  not  containing  zv,  there 
is  a  round  cut  RjjU  satisfying  |A(i?)|  <  £i|A(£/)|. 

Proof.  Set  =  Ki  +  b~1,/2a!-1  and  £1  =  4A f[(/3  +  e~l).  Let  U  C  G'v  be  a  set  not  containing 
zv,  and  write  p  =  | <5(C/) | ;  we  construct  a  round  cut  R  containing  U  for  which  |A(i?)|  <  fip. 
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Now  if  we  contract  G'V\U  to  a  single  vertex,  we  obtain  a  planar  graph  with  maximum 
face  size  K\  (as  opposed  to  £;  this  is  due  to  the  large  spacing  of  the  vertices  of  a'v ).  So  by 
Lemma  4.7,  the  maximum  distance  between  two  points  on  tt(U)  is  at  most  K\p. 

Next  we  claim  that  every  vertex  in  U  must  be  within  distance  a~1/2a~1p  of  n (U).  The 
reason  for  this  is  analogous  to  the  proof  of  Lemma  4.6  —  if  not,  then  U  would  contain  a 
ball  of  more  than  this  radius,  which  would  contain  more  than  a-1p2  vertices;  a  contradiction 
since  |h(I/)|  =  p. 

Thus  for  any  u  G  7r(I/),  we  have  U  C  Br(u),  where 

r  =  (/«i  +  a^l^2a^l)p  =  £[■ p . 

Now  by  Lemma  4.5,  there  is  an  r'  between  r  and  2 r  such  that 

l<K-Br'(u))|  <  4/3  A f[p. 

Now  let  R  D  U  denote  the  round  cut  Br>(u )  ft  Cv.  Every  edge  of  S(R)  is  an  edge  of  either 
5(Bri(u))  or  of  S(CV  fl  R ).  The  former  quantity  was  just  shown  to  be  at  most  4/3A£jp.  To 
bound  the  latter  quantity,  note  that  any  two  vertices  in  n (Cv  fl  R)  are  at  most  2 r'  <  4 C,[p 
apart;  since  the  facial  cycle  containing  n(Cv)  is  e-smooth,  this  means  that  n(Cvr\R)  contains 
at  most  4e-1£jp  vertices,  and  hence 

^(ani?)!  <  4Ae_1^p. 

The  claim  now  follows  since 

\S(R)\  <4Aei(/3  +  e->  =  eiP- 


Thus  we  have 

Theorem  5.18  There  is  a  randomized  (off-line)  mdpp  algorithm  in  uniformly  densely  em¬ 
bedded  Eulerian  graphs  that  produces  a  constant-factor  approximation  with  high  probability. 


6  Extensions 

6.1  Durations  and  Profits 

In  the  on-line  algorithm  we  assume  that  (i)  all  connections  have  infinite  duration,  and  (ii)  all 
connections  have  the  same  “value”  (i.e.  our  objective  function  could  have  been  a  weighted 
sum  of  the  set  of  pairs  we  accept,  rather  than  an  unweighted  sum).  However,  there  are 
general  transformation  techniques  due  to  Awerbuch  et.  al.  [5]  that  allow  us  to  convert  our 
results  to  on-line  algorithms  that  can  handle  connections  of  limited  duration  and  variable 
value,  at  the  cost  of  additional  logarithmic  terms  in  the  competitive  ratio.  Specifically,  we 
pay  0(log T)  and  O(logP),  where  T  and  P  are  the  ratios  between  the  largest  and  smallest 
durations  and  values  respectively. 
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6.2  Graphs  with  an  Exceptional  Face 

In  this  section,  we  sketch  the  extension  of  our  algorithms  to  densely  embedded,  nearly- 
Eulerian  graphs.  Recall  from  Definition  4.3  that  such  a  graph  satisfies  the  properties  of  a 
uniformly  densely  embedded  Eulerian  graph,  except  that  it  is  allowed  to  contain  an  “excep¬ 
tional”  face  <&*,  with  facial  cycle  Q*  that  may  have  length  greater  than  £  and  may  contain 
vertices  not  of  even  degree. 

For  the  remainder  of  this  section,  let  G  denote  a  densely  embedded,  nearly-Eulerian 
graph  with  parameters  a,  A,  A,  and  £.  For  simplicity,  we  assume  that  the  facial  cycle  Q*  is 
sufficiently  large  that  it  is  not  contained  in  any  set  B\\ogn(v):  it  is  not  difficult  to  remove 
this  assumption. 

The  changes  required  in  the  algorithm  come  from  the  fact  that  there  can  now  be  a  G- 
normal  curve  joining  two  distant  vertices  in  G  that  intersects  G  relatively  few  times  —  this 
is  because  it  can  pass  through  the  large  face  <&*.  This  has  consequences  in  the  proofs  of 
Lemma  5.6  (and  its  relatives)  and  Claim  5.13.  However,  by  requiring  the  outer  cycles  of 
clusters  and  enclosures  to  satisfy  a  more  restrictive  notion  of  e-smoothness,  these  facts  will 
follow  as  before. 

We  define  our  more  restrictive  type  of  e-smoothness  as  follows.  Let  G/Q*  denote  the 
graph  G  with  a  single  additional  node  q*  joined  by  length-0  edges  to  each  vertex  of  the  long 
facial  cycle  Q*.  Then  a  small  cut  passing  through  two  distant  vertices,  as  described  in  the 
previous  paragraph,  does  correspond  to  a  short  path  in  G/Q*  —  it  simply  makes  use  of  the 
additional  node  q*.  Now  it  is  straightforward  to  show  that  we  need  only  require  the  outer 
cycles  of  the  clusters  and  enclosures  to  be  e-smooth  in  the  graph  G/Q*]  and  this  can  be 
accomplished  by  running  the  e-smoothing  algorithm  in  this  graph  instead  of  in  G. 

This  introduces  a  further  difficulty,  however.  Say  that  we  have  just  smoothed  some  cycle 
Q ,  obtaining  a  cycle  Q' .  While  Q'  will  be  e-close  to  the  original  cycle  Q  in  G/Q*,  there  is 
no  reason  why  this  means  it  will  be  e-close  in  G. 

To  handle  this,  we  strengthen  the  statement  of  Theorem  5.4,  as  follows.  For  a  vertex  u, 
define  the  u-restricted  distance  du(v,w )  between  two  vertices  v  and  w  to  be  the  minimum 
length  of  a  v-w  path  avoiding  u.  From  the  proof  of  Theorem  5.4,  one  sees  that  we  can 
always  find  a  short  path  from  the  smooth  cycle  Q'  back  to  the  original  cycle  Q  that  avoids 
any  prescribed  vertex;  i.e.  Q'  is  e-close  to  Q  with  respect  to  any  u-restricted  distance  function. 
In  particular,  Q'  is  e-close  to  Q  in  G/Q*  with  respect  to  the  (^-restricted  distance  function; 
that  is,  Q'  is  e-close  to  Q  in  the  original  graph  G. 

Thus,  we  can  obtain  clusters  and  enclosures  with  boundaries  which  are  e-smooth  in  G/Q*, 
and  which  are  not  far  from  the  original  boundaries  in  G.  The  sets  t'v  will  now  consist  of 
evenly  spaced  vertices  only  on  the  part  of  an  enclosure’s  outer  facial  cycle  that  does  not  lie 
on  <f>*.  The  proof  of  Lemma  5.6  now  follows  exactly  as  before.  When  we  use  Schrijver’s 
theorem  to  construct  inter-cluster  paths,  we  also  cut  the  surface  E  along  the  long  facial  cycle 
Q*  so  as  to  remove  from  the  surface.  Now  an  essential  curve  can  be  anchored  on  the 
boundary  of  as  well  as  on  the  boundary  of  an  enclosure;  but  this  poses  no  problem  since 
the  enclosure  boundaries  are  e-smooth  in  G/Q*. 

Once  the  simulated  network  J\f  has  been  set  up,  the  on-line  and  off-line  algorithms  work 
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exactly  as  before.  (In  particular,  Lemma  5.17  follows  without  modification,  since  the  large 
face  $*  is  incorporated  into  the  hypotheses  of  Lemma  4.7.)  Thus  we  have, 

Theorem  6.1  There  is  an  O (log n)- competitive  on-line  mdpp  approximation  in  any  densely 
embedded  nearly-Eulerian  graph. 

Theorem  6.2  There  is  a  randomized  (off-line)  mdpp  algorithm  in  any  densely  embedded 
nearly-Eulerian  graph  that  produces  a  constant-factor  approximation  with  high  probability. 
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